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ABSTRACT 


Combat  models  often  involve  target  detection  times  which  may  vary  with 
different  observers  due  to  characteristics  of  personnel,  or  detection  systems.  They 
may  also  be  affected  by  different  environmental  factors  such  as  visual  levels,  sea 
states,  terrains,  etc.  There  is  often  interest  in  quantifying  the  effects  of  different 
observer  characteristics  and  environmental  factors  on  detection  times.  A  hierarchical 
gamma/Weibull  regression  model  is  considered  which  can  incorporate  observer 
characteristics  and  environmental  effects  which  may  influence  the  time  to  detect 
targets.  Numerical  procedures  for  the  estimation  of  parameters  of  the  hierarchical 
gamma/Weibull  model  based  on  maximum  likelihood  are  described.  Results  of 
simulation  experiments  to  study  small  sample  behavior  of  the  estimates  are  reported. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Combat  models  often  involve  target  detection  times  which  may  vary  with 
different  observers  due  to  characteristics  of  personnel,  or  detection  systems.  They 
may  also  be  affected  by  different  environmental  factors  such  as  visual  levels,  sea 
states,  terrains,  etc.  There  is  often  interest  in  quantifying  the  effects  of  different 
observer  characteristics  and  environmental  factors  on  detection  times.  In  Chapter 
III  a  hierarchical  gamma/Weibull  regression  model  [Ref.  1]  is  considered  which  can 
incorporate  observer  characteristics  and  environmental  effects  which  influence  the 
time  to  detect  targets.  Numerical  procedures  for  the  estimation  of  parameters  of  the 
hierarchical  gamma/Weibull  model  based  on  maximum  likelihood  are  described. 
Results  of  simulation  experiments  to  study  small  sample  behavior  of  the  estimators 
are  reported.  In  the  remainder  of  this  chapter,  two  models  related  to  the  one 
considered  in  this  thesis  will  be  described.  Numerical  procedures  for  parametric 
estimation  for  these  models  are  studied  in  Bae  [Ref.  2].  Experience  gained 
concerning  the  numerical  stability  of  these  procedures  will  be  described. 

B.  TWO  RELATED  MODELS 

Two  parametric  models  for  the  distribution  of  time  to  detect  targets  were 
considered  in  a  thesis  written  by  a  Korean  officer, Mr.Bae  in  Sept., 1989  [Ref.  2]. 

Model  1  :[Ref.2:pp.2-14]  supposes  that  there  are  M  observers.  Observer  i  has 
parameter  A,  which  characterizes  the  ability  of  the  observer  to  detect  a  target. 
Observer  i  is  presented  with  N;  targets.  The  parameters  of  all  M  observers, 
Aj,A2,...,Am,  are  assumed  independent  with  a  common  two  parameter  gamma 
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distribution.  Given  A;,  the  detection  times  of  observer  i  are  conditionally 
independent,  having  Weibull  distributions  with  known  shape  parameter  e'*1  and 
known  scale  parameter  for  j  =  l,2,..,Ni.  Some  of  the  detection  times  may  be 
censored.  Bae  [Ref.  2:pp.  2-20]  reports  on  numerical  procedures  to  find  the 
maximum  likelihood  estimates  for  the  shape  and  scale  parameters  of  the  gamma 
distribution. 

Model  2  :  is  a  Weibull  regression  model  [Ref.2:pp.21-30],  It  assumes  observer 

i  has  explanatory  variables  x^,  xij2, . .  xijp  relating  to  his  j,h  target  representing 

factors  which  influence  his  time  to  detection.  The  detection  times  for  the  observers 

are  independent  random  variables  having  Weibull  distributions.  Again,  some  of  the 

detection  times  may  be  censored.  The  scale  parameter  of  the  Weibull  distribution  for 

the  detection  times  of  the  j,h  target  by  observer  i  is  of  the  form  where 

p  , 

i  % P*  ;  the  shape  parameter  is  of  the  form  e  1  .  Bae  [Ref.  2:pp.  21-36] 

*  *»i 

reports  on  numerical  procedures  to  find  the  maximum  likelihood  estimates  of  {ft} 
and  {ft.}. 


C.  OBJECTIVE  AND  METHODOLOGY  IN  THE  COMBINED  MODEL 

The  main  effort  of  this  thesis  will  be  to  study  a  full  hierarchical  gamma/ Weibull 
regression  model  which  is  a  combination  of  the  two  previous  models.  Since  the 
models  considered  in  Bae  [Ref.  2]  will  be  special  cases  of  the  one  considered  in  this 
thesis,  it  is  expected  that  instabilities  found  in  the  numerical  procedures  to  estimate 
the  parameters  in  the  two  previous  models  will  appear  for  the  combined  model.  Two 
numerical  instabilities  that  result  from  the  numerical  procedures  of  Bae  [Ref.  2]  are 
described  below. 

In  Model  1,  a  modified  Newton-Raphson  method  [Ref.  3]  is  applied  to  solve 
the  nonlinear  system  of  equations  for  the  maximum  likelihood  estimates  of  the 
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gamma  parameters.  It  is  found  that,  particularly  for  small  sample  sizes,  the  slope  of 
the  likelihood  surface  is.  very  flat  resulting  in  numerical  instability  for  estimation  of 
the  scale  parameter  i).  As  a  result,  in  the  full  hierarchical  model  presented  here  we 
have  used  a  single  parameter  gamma  distribution  as  the  second  stage  model.  The 
maximum  likelihood  estimate  of  the  parameter  is  found  by  the  bisection  root  search 
method  [Ref.  4].  • 

A  numerical  overflow  phenomenon  was  occasionally  found  in  the  Weibull 
regression  model  while  iterating  the  Newton  procedure  for  {£;}  in  Model  2 
[Ref.2:p.25].  Once  again,  the  flatness  of  the  relationship  determining  the  estimate 
of  {£;}  resulted  in  this  instability.  As  a  result,  the  Newton  procedure  to  determine 
{£;}  in  the  gamma/Weibull  regression  model  of  this  thesis  has  some  checks  to  detect 
this  numerically  instability. 

In  Chapter  II,  a  numerical  procedure  is  presented  to  find  the  maximum 
likelihood  estimate  for  a  single  parameter  gamma  version  of  Model  1.  In  Chapter  III, 
the  full  hierarchical  model  is  presented  and  numerical  procedures  to  estimate  the 
parameters  given.  In  both  chapters,  simulation  results  are  presented  to  study  the 
small  sample  behavior  of  the  estimates.  All  simulations  were  carried  out  on  an  IBM 
3179  G  mainframe  computer  at  the  Naval  Postgraduate  School  using  the  APL 
GRAFSTAT  random  number  package  [Ref.  5], [Ref.  6]. 

We  hope  that  the  hierarchical  gamma/Weibull  regression  model  and  these 
estimation  procedures  will  be  a  useful  tool  to  describe  and  predict  target  detection 
times,  which  are  one  aspect  of  the  effect  of  human  performance  on  the  battlefield. 
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II.  A  SINGLE  PARAMETER  GAMMA  HIERARCHICAL  MODEL 
A.  MODEL  DESCRIPTION 

This  model  is  very  similar  to  Model  1  in.  [Ref.  2]1.  The  difference  is  that  the 
parameters  Aj,  i  =  which  reflect  the  abilities  of  observer  i  to  detect  targets, 

are  assumed  to  be  independent  having  a  single-parameter  gamma  distribution \, 
GAM( a,a)  rather  than  having  two  parameters.  The  form  of  gamma  density  function 
used  here  is 

$(0)  =  f0r  6>o.  (2.1) 


For  numerical  reasons  we  will  parametrize  a  as  e\  For  convenience  we  will 
reiterate  the  remainder  of  the  assumptions  of  Model  1  given  in  [Ref.  2];  Given  Aj 

=  9  ,  the  times  it  takes  for  observer  i  to  detect  target  j,  j  =  1,2 . Ni(  denoted  by  Ujj, 

are  assumed  to  be  conditionally  independent  random  variables  with  Weibull 
distributions ,  WEli^^e^1^  having  cumulative  distribution  function 

-  expf-ef#/'1}  no  (2-2) 

independent  of  other  observers.  The  variations  of  the  A,  are  introduced  to  represent 
the  individual  differences  between  the  observers.  When  the  ith  observer  is  presented 
’  with  his  jth  target,  only  an  opportunity  time  Ojj  is  allowed  for  him  to  detect  it.  An 
observer  either  successfully  detects  the  target  within  this  time  or  never  detects  it. 

1  In  the  model  of  [Ref.  2],  a  two  parameter  gamma  distribution  is  used. 
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Data  for  the  ith  observer  consist  of  times  of  detection  for  the  successes  and  the 
lengths  of  opportunity  times  for  the  failures.  For  each  i  =  1,2,...,M  ,  j  =  1,2, ...,N,  , 
let 

Yi}  =  min(lnl/.,lnO..)  (2.3) 


and 


1  ifU^Oy 
0  otherwise 


(2.4) 


The  Yy  are  the  censored  ln-detection  times  and  is  an  indicator  of  whether  the  In¬ 
time  to  detect  the  jlh  target  by  the  ith  observer  is  censored  or  not.  Let 


N, 

C.  =  E  A.. 

i  “i ij 


(2.5) 


be  the  number  of  targets  detected  by  observer  i.  In  this  chapter,  we  will  assume  the 
Weibull  parameters  {/i^}  and  {£,}  are  known  constants,  and  we  are  only  interested 
in  estimating  the  gamma  parameter  t].  In  the  next  section,  the  likelihood  function 
for  rj  is  given  and  the  bisection  method  for  finding  the  maximum  likelihood  estimate 
of  rj  is  described.  In  the  final  section,  we  present  results  of  a  simulation  study  of  the 
behavior  of  the  estimators  for  small  sample  sizes. 
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B.  MAXIMUM  LIKELIHOOD  ESTIMATION  AND  THE  BISECTION  ROOT 
SEARCH 


1.  The  Likelihood  Function 

Given  Aj  =  6,  the  conditional  likelihood  function  [Ref.2:p.3j  for  observer 
fusing  the  censored’  ln-times  Yjj  is 


y-i 


(2.6) 


Let 


S : 


/■I 


(2.7) 


and 


n, 


K, ■  =  expj  EA|yrln^)e'£i 


(2.8) 


Equation  (2.6)  can  be  rewritten  as  follows : 

i, <ii, lie) » e^expj-es,) .  (2-9) 


From  the  equations  (2.1)  and  (2.9),  it  follows  that  the  unconditional  likelihood 
function  for  observer  i  is 
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(2.10) 


*  /c'£/ii.i!0is(exio 

.  f-eV«,a(«9)-1g-V9 

J0  f(«) 


=  £ 


a 


a  Crl 


‘  (Si+u)Ci+a  k-0 


n  (o+*) . 


C,-l 

In  equation  (2.10),  if  Q  =  0,  then  ]~[  (a+k)  =  1  .  Recall  the  parametrization  a 

*»o 

=  e”.  The  unconditional  ln-likelihood  function  for  observer  i  can  be  rewritten  as 


lnZ/li&Ty)  = 

In#,  +  -  (C<+e’>)lh(Sj+e,>) 


c,-i 

+  Eln(en+fc) 

i=0 


(2.11) 


cri 

where  if  Q  =  0,  then  £  ln(e ’>+/:)  =  0  .  Since  the  observers  are  independent,  the 

*=o 

unconditional  ln-likelihood  for  all  M  observers  is 


M 

In L  =  lnl,  = 

*  ,  w  Crl  1  (2.12) 

T  jln/fi+enr|-(C|.+en)ln(S(.+e,1)+  E  ln(e,i+*)|. 
i=i  l  *=o  J 


In  this  chapter,  we  will  focus  attention  only  on  finding  the  maximum  likelihood 
estimate  of  r\.  The  derivative  of  the  ln-likelihood  with  respect  to  t]  is 
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ainL 

Sti 


■  /Oi) 


M 

M 


Cr  1 


S-C: 

r\  +  -  In (S.+eT1)  +  E 


1 


S.+eT> 


(2.13) 


The  maximum  likelihood  estimate  for  r\  is  the'-solution  of  the  equation  f(rj)  =  0. 

2.  Bisection  Root  Search  for  r? 

Initially,  the  Newton  procedure  was  used  to  numerically  solve  the  equation 
f (77)  =  0.  Unfortunately,  the  procedure  frequently  either  converged  to  an 
unreasonable  number  or  encountered  problems  of  numerical  instability.  Plots  of  f(rj) 
indicate  that  f  can  have  multiple  zeros,  ^representative  plot  of  f  is  shown  in  Figure 
A  2  (next  page).  In  Figure  A,  one  v  symbol  indicates  a  zero  of  f,  while  two  v 
symbols  indicate  the  reasonable  root  to  be  used  as  the  estimate  for  r?. 


•y 

“  This  is  a  typical  graph  for  f(»j)  which  was  generated  for  simulated  data  with  0  =  10,  M = 15,  N  =  15 
and  random  seed  =  16807.  The  true  value  of  1 7  is  1. 
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Figure  A  Multiple  Roots  of  the  Derivative  of  the  In-Likelihood  Function  for  rj 


As  a  result  of  these  possible  multiple  zeros,  the  bisection  root  search  method  [Ref.4] 
is  used  to  solve  f(rj)=0.  The  following  sub-section  will  detail  the  estimation 
procedure  for  tj. 
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;  a.  Initial-Estimate  for  r? 

The  initial  estimate  for  rj  is-similar  to  its  initial  estimate  in  the  two- 
parameter  gamma  model  [Ref.  2],  and  is  obtained  by  setting  3  /30  =  0  in  that 
procedure.  Consequentially,  the  initial  estimate  of  r\  is 


i°  =  -lnl 


-x  »2  ir 


(2.14) 


where 


m,  = 


= 


-0.5772 

El 

l»V«l 


M  N, 

EEA(. 

i-V«l  ‘ 


(2.15) 


If  <  0  ,  set  tj°  =  0 4. 
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b.  Determination  of  Bounds  for  the  Root 

The  purpose  of  this  step  is  to  bound  the  root  that  will  be  used  as  the 
estimate  between  two  points,  called  t?L  and  t?n.  One  of  these  points  is  t?°  and  the 
other  is  found  as  follows : 


3  In  MriBae’s  model  [Re£2],  the  regular  gamma  distribution  is  parametrized  with  shape  parameter  v  =en 
and  scale  parameter  a=e\  Setting  /J0=0,  a  =  u  =  cn  yields  the  one-parameter  gamma  distribution. 


4  The  estimate 


A  2  A  2  7C  .  A 

nu-mx - <0 

2  1  6 


suggests  that  the  variability  of  the  data  is  small.  Hence,  it 


might  be  better  to  set  rf  equal  to  a  large  positive  number,  such  as  13.  A  topic  for  future  research  is 
to  explore  better  initial  values  for  rj  in  this  default  case. 
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r\L=T\°-kx5 

Set  .  Oft  „  where  5  is  a  small  increment. 

T)/j=1i  +X'xo  je=ly2,3 . 

Compare  the  signs  of  f(tjL)  and  f (t?r)  with  f(7?°)  until  the  first  different  sign  from  f(r?0) 
is  . found.  This  yields  a  bound  on  the  closest  root  to  ??°. 
c.  Bisection  Search 

After  determining  the  left  bound,  r?L  and  right  bound  ,r/R,  we  begin 
the  bisection  search  by  picking  the  midpoint  between  them,  called  ?]M,  and  computing 
the  value  of  f (r?M).  If  f(rjM)  and  f(rjL)  have  the  same  signs,  the  root  is  on  the  other 
side  of  r?M;  thus,  move  the  left  bound  r?L  to  midpoint  r?M.  If  the  signs  are  different, 
move  the  right  bound  to  the  midpoint.  Iterate  the  bisection  search  until 
- <0.01  .  The  maximum  likelihood  estimate  of  r?  is  r?M, 

C.  IMPLEMENTATION  OF  THE  SIMULATION  AND  RESULTS 
1.  Simulation 

In  order  to  study  the  sampling  properties  of  the  maximum  likelihood 
estimate  of  r?,  different  numbers  of  observers  (M),  targets  (Nj),  and  opportunity  times 
(Ojj)  are  used  as  inputs  to  the  simulation.  Each  set  of  inputs  (M.Nj.Ojj)  has  a 
simulation  run  of  (R)  100  replications  and  each  run  begins  with  the  same  initial 
random  seed  (466801743).  Statistics,  of  mean  bias  (M.B),  mean  square  error  (M.S.E) 
and  their  standard  errors  (  S.E(M.B)  and  S.E(M.S.E)  ),  are  computed  according  to 
following  definitions : 

=  4EM,-n)  (2-16a) 
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Ct 


s:e(M.b )  =  ; 

^ 


m.s,e  =  |E(nr^)2  (2-17a) 

/tr«i 


S.E(MXE)  = 


(2.17b) 


1 


R(R-Z)r 


r-x\)-M.Bf 


(2.16b) 


where  f[  r  represents  the  point  estimate  of  parameter  r?  in  the  r,h  replication.  The 
fraction  of  detection  level  or  the  averaged  uncensoring  level  (UC)  for  R  (100) 
replications  is  determined  by 


UC 


R  M 

EE 

r=l  hi 
R  M 

E  EW 

r«l  i»l 


(2.18) 


where  C;(r)  is  the  number  of  targets  detected  by  observer  i  in  the  r,h  replication  and 
Nj(r)  is  the  number  of  targets  presented  to  observer  i  in  the  r,h  replication. 
Theoretically,  the  longer  the  opportunity  times,  the  higher  the  UC  level  should  be. 
An  outline  of  a  replication  in  the  simulation  is  as  follows : 
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a.  Input  to  Simulation 


•  M  the  number  of.  observers  (  Possible  values  are  15 , 25  and  35  ) 

•  Nj  the  number  of  targets  presented  to  the  ith  observer ,  where  i  = 

1,2....M  (  Possible  values  are  15  ,25  and  35) 

•  Ojj  the  opportunity  time  for  ith  observer  to  detect  the  jth 

target, (Possible  values  are  10  and  15).  All  of  the  opportunity 
times  are  the  same. 

•  (j  the  shape  parameter  of  the  Weibull  distribution  for  ith  observer, 

where  i  =  1,2, ...M  ;  (  Set  =  0  for  all  i ) 

•  /ijj  the  scale  parameter  of  the  Weibull  distribution  ,WEI(/tn ,  e'5i); 

(Set  /ijj  =  4.2  ) 

•  T)  the  parameter  for  the  single-parameter  gamma  distribution  , 

GAM(tt,a) ,  with  a  =  en  ;  (Set  the  true  value  of  rj  =  1 ) 


b.  Simulation  of  Data  for  ith  Observer 


Generate  independent  single-parameter  gamma  r.v.’s  Aj  from 
GAM(at, a), having  the  density  function  as  in  equation  (2.1). 


Generate  independent  exponential  distribution  r.v.’s  Wi;  with  mean 
1  i.e.  Wjj  ~  EXP(l)  for  all  i  and  j  J 

( W,  V 

Compute  the  target  detection  times :  U{j  =  p J 


.V  J 

Compute  the  recorded  ln-times  Y;:  as  in  equation  (2.3). 


Compute  Aj:  as  in  equation  (2.4). 


’ 

(WB) 

«'<  ] 

f 

( t  \ 

V 

s/  M,-e 

►  =  pIw^q 

conditional  Weibull  distribution  of  equation  (2.2). 


1-exp' 


/  \ 

-U 

e 

i 

B 

[hi 

which  is  the 
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c.  Initial  n  and  Bisection  Search 


•  Compute  m2  as  in  equation  (2.15)  and  determine  the  initial  value 
of  rj  as  in  equation  (2.14). 

•  Perform  the  Bisection  root  search  procedure  for  the  maximum 
likelihood  estimate,  of  rj,  shown  as  the  flowchart  in  Figure  B. 


All  these  procedures  were  written  in  APL  codes, named  "SIMULA1", 
"GAMMA",  and  "FVALUE",  which  are  listed  in  Appendix  C. 


14 


Aljroritltm  of  Bisection  Search 


Figure  B  The  Algorithm  of  Bisection  Search 
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2.  Results 

The  results  of  the  simulation  experiments  appear  in  Tables  and  graphs  in 
Appendix  A.  Tables  1  and  2  contain  the  general  statistics  for  the  bias  of  p  at  the 
opportunity  times  O  =  10  and  15  respectively.  Table  3  in  combination  with  box 
plots,  shown  in  Figures  1  through  3,  presents  statistics  of  the  bias  of  the  estimates  of 
T)  as  a  function  of  the  numbers  of  targets  and  observers  in  the  case  O  =  10.  Figures 
4  to  6  show  histograms  of  the  bias  with  the  number  of  targets  fixed  and  the  number 
of  observers  changing.  Similarly,  Table  4  and  Figures  7  through  12  focus  on  the  case 
O  =  15.  In  summary,  the  simulations  indicate  the  following  : 


•  Opportunity  time  (0) 

Tables  1  and  2  indicate  that  a  longer  opportunity  time  is  reflected  in 
a  higher  uncensoring  level  (UC),  which  slightly  decreases  the  mean 
bias  (M.B)  and  mean  square  error  (M.S.E)  of  r\. 

•  The  Number  of  Observers  (M) 

Increasing  the  number  of  observers  for  a  fixed  number  of  targets 
results  in  greater  changes  in  the  mean  bias  and  mean  square  error  of 
T]  than  increasing  the  number  of  targets  with  a  fixed  number  of 
observers.  All  the  box  plots  and  histograms  display  this  tendency. 
This  phenomenon  is  quite  reasonable,  because  the  parameter  r?  in 
gamma  distribution  reflects  the  variability  in  abilities  of  observers. 
Hence  data  with  more  observers  will  provide  better  estimation  of  r? 
( less  bias  and  standard  error ). 

•  The  Number  of  Targets  (Nj) 

The  box  plots  indicate  that  increasing  the  number  of  targets  for  a 
fixed  number  of  observers  has  a  small  effect  on  M.B  and  M.S.E  of  r? 
with  no  systematic  trend. 
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lit  THE  HIERARCHICAL  GAMMA /WEI BULL  REGRESSION  MODEL 
A.  MODEL  DESCRIPTION 

Recalling  the  assumptions  of  the  model,  we  suppose  there  are  M  observers, 
indexed  by  i.  The  ith  observer  is  presented  with  Nj  targets,  indexed  by  j.  The 

variable  is  the  time  it  takes  for  observer  i  to  detect  the  target  j.  Let  x^  ,  Xjj2 . . 

Xjjp  be  the  values  of  environmental  variables  (e.g.  terrain,  sea  state  ,  atmospheric 
condition  etc.)  which  may  affect  U^.  Given  the  quantified  value  of  the  ability  for 
observer  i,  Aj  =  6,  the  Ujj’s  are  assumed  conditionally  independent  random  variables 
having  Weibull  distributions,  ,  where  the  scale  parameter,  instead  of 

being  constant,  has  the  form 

Vij  =  «*{*<,£■}  =  exp{^P*}-  (3,1) 


The  regression  coefficient  is  the  contribution  to  the  Weibull  scale  parameter  of 
the  k,h  explanatory  variable.  /50  is  usually  the  constant  term  in  the  regression.  The 
random  variables  {A,}  are  assumed  independent  from  a  single  parameter  gamma 
distribution,  GAM  (a,  a)  with  a  =  en. 

B.  MAXIMUM  LIKELIHOOD  ESTIMATES  FOR  ALL  PARAMETERS 

The  model  assumptions  detailed  above  are  the  same  as  those  in  Chapter  II, 
except  that  the  scale  parameter  gl}  in  the  Weibull  distribution  is  now  a  function  of 


17 


the  explanatory,  variables.  In  this  section  we  consider  the  estimation  of  all  the 
parameters  of  r? ,  (4)  and  {^}.6 

From  equations  (2.7),  (2,8),  (2.12)  and  equation  (3.1),  we  can  rewrite  the 
unconditional  ln-likelihood  function  for  all  M  observers  as  follows : 


InL  * 

m  [  -  cr 1 

g{r)£&)  «  £  +  -  (C,.+^)ln(S.+e»)  +  Eln(eVJt) 

i-l  l  - 


Yjj ,  Ay ,  and  Q  are  defined  in  equations  (2.3),  (2.4)  and  (2.5)  respectively. 

In  the  following  sub-sections,  we  will  describe  the  procedures  which  are 
employed  to  estimate  all  the  model  parameters.  The  general  procedure  includes  the 
initial  rough  estimation  of  (4)  and  {£j}  without  the  hierarchical  gamma  r.v.’s;  a 
variability  check  for  the  need  to  include  the  hierarchical  gamma  parameter  r?;  and 
finally,  a  large  recursive  procedure  for  estimating  all  the  model  parameters  if  a 
hierarchical  model  is  necessary. 

1.  Initial  Estimation  of  a  Simple  Weibull  Regression  Model 

An  initial  model  for  the  data  is  the  simple  Weibull  regression  model  for 
the  target  detection  times  described  previously  as  the  Model  2  of  Bae  [Ref.  2:pp.  21- 
30].  The  key  steps  of  the  estimation  procedure  are  summarized  as  follows : 

6  In  the  Model  2  of  Bae  [Ref.  2],  no  gamma  variable  was  involved.  Only  the  Weibull  parameters 
{&}  and  {£d  were  considered. 
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Step-1.1  :  Initially,  set 7  2Ljl0  =  In Utj  and  =£•  =  0  . 

Step-1.2  :  This  is  an  iteratively  re-weighted  regression  step;  cf.  McCullagh 
and  Nelder  [Ref.  7]. 


u  -  [y«  -  */)< 


WU  =  \/exP(r</}  » 


(3.3) 


(3.4) 


=  XiJh'Wije 


-V 


where  h=0,l,..,p  (3.5) 
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(3.6) 


Regress  the  dependent  variable  Zjj  against  the  independent  variables  uijh, 
for  h  =  0,1,..., p  ;  e.g.  the  estimate  has  the  form 

B  =  (3-7) 


where 


7 


From  this  page  on,  all  parameters  with  superscript  0  imply  the  current  values  of  their  estimates. 
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(3.8), (3.9) 


Step- 1.3  :  Compute  the  new  £,  by  an  approximation  based  on  Newton’s 
procedure.  First,  update  the  value  of 


r,j  =  {ytJ  ~  P£..P.°)e'('  ,  where  &°  =  U  . 


(3.10) 


Compute 
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Compute-  the  new  value  of  2  B°  =  E  p°  and  iterate  the  regression, 

"  k*  0 

starting  at  Step- 1.2.  for  one  more  time. 

2.  Does  the  Simple  Weibull  Regression  Adequately  Account  for  the 
Variability  of  the  Data  ? 

The  full  hierarchical  model  will  be  used  if  the  data  indicate  more 
variability  than  can  be  described  by  the  simple  Weibull  regression.  If  the  data  do 
come  from  the  simple  Weibull  model,  then 
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will  have  the  same  distribution  as  independent  exponential  random  variables8  with 
mean  1.  To  assess  the  variability  not  accounted  for  by  the  simple  regression  model, 
we' compute 


A  2 

°E  = 


and 


—L-V 

(N-l)r 


-1 


M 

N='£Ni  . 

i- 1 


(3.12) 


♦  1 

Since  the  variance  of  a  unit  exponential  is  1,  if  d  E  <  1  ,then  we  will  stop  the  model 
fitting  with  the  simple  Weibull  regression  model.  If  a\  >  1,  we  go  on  to  the 
hierarchical  model  to  explarn  the  extra-variability.  Unfortunately,  crude  estimate  of  b\ 
uses  both  the  observed  and  unobserved  Uy.  A  topic  for  future  research  is  to  improve 
the  estimate  of  the  variability  not  accounted  for  by  the  simple  regression  model. 

Step-1. 5  :  If  bE  $  1 ,  stop  and  use  the  Weibull  regression  model  only. 
Otherwise,  the  data  will  be  modeled  with  the  full  hierarchical  model.  A  recursive 
procedure  to  estimate  t?  ,  {/)k}  and  {^}  follows  and  is  initialized  with  and 

computed  above. 


8  .)''''}= 1  - e' W  which  is  the  pure  Weibull 

distribution. 
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3.  A  Recursive  Estimation  Procedure  for  the  Parameters  of  the  Full 
Hierarchical  Model 

,  a.  Bisection  Search  with  Instability  Check  for  Estimating  r? 

The  partial  derivative  of  equation  (3.2)  with  respect  to  ??>  -  ,  is 

3ri 

the  same  as  the  derivative  of  the  ln-likelihood  with  respect  to  r)  for  the  model  of 
Chapter  II.  The  derivative  appears  in  (2.13)  and  once  again  the  problem  is  to  solve 
the  equation  f(r?)=0.  Using  the  estimates  and  an  initial. estimate  for 
r?  is  computed  as  in  equation  (2.14) ;  that  is , 


=  -Info 


.2  .2  n 
m2-mr— 


(3.13) 


where 


(3.14) 


If  <0  ,  set  r?°  =  0. 9 

6 

In  addition  to  determining  the  bounds  for  a  reasonable  root  to  be  the 
estimate  of  r\  as  described  in  Chapter  II,  a  check  for  the  shape  of  f (t?)  is  also 
performed.  This  additional  check  is  performed  because  two  of  the  simulation 
replications  exhibited  a  very  flat  f(t?).  The  graph  of  f(rj)  for  one  of  the  two 


9  Same  as  footnote  4  on  page  10. 
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replications  appears  in  Figure  C  .  This  unusual  shape  for  f (rj)  results  in  the 
bisection  method  converging  to  an  unreasonable  foot.(e:g.  r\  =146  in.the  Figure  C). 
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Figure  C  An  Ill-Conditioned  Derivative  of  the  ln-Likelihood  Function  for  r? 

The  graph  in  Figure  C  indicates  that  the  values  of  f(rj)  are  asymptotic 
to  a  fixed  nonzero  value  as  rj  tends  to  -®;  this  behavior  is  the  same  as  in  Figure  A. 
However,  when  rj  is  positive,  f  has  very  flat  curve.  The  values  of  f(rj)  are  positive 


u  The  graph  is  one  of  the  only  two  ill-conditioned  cases  for  f(rj)  in  the  simulations.  This  case 
occurred  in  the  ganuna/WeibuIl  regression  model  with  0=10,  M  =  15,  N=25  and  random  seed  = 
1905665785.  The  other  case  occurred  with  0  =  15,  M  =  15,  N=25  and  same  random  seed.  The  true  value 


and  close  to  zero  but  the  curve  does  not  actually  cross  zero  until  rj  is  large.  The 
following  check  is  done  to  detect  this  situation.  If  a  reasonable  right-hand  bound  for 
a  zero  of  f (rj)  cannot  be  found,  check  the  values  of  f(r?)  and  f '(??)  for  r?  positive  if 
the  absolute  value  of  f(r?)  is  less  the  ex  and -the  absolute  value  of  f'(r?)  is  lesa  than 
e2,  where  1  >  >  >  e2  >  0,  then  the  simulation  stops  for  -this  set  of  random 

numbers  arid  new  random  numbers  are  generated.  Otherwise,  perform  the  bisection 
search  as  described  in  Chapter  II.  In  the  procedure,  f'(r?)  for  rj  positive  is  estimated 
as 


where  **»  =  ’1"tto<6 

nw.1)=<l0H*-l)x6. 


(3.15) 


The  algorithm  of  this  bisection  search  with  the  check  for  flatness  of  the  curve  of  f 
will  be  detailed  later  in  Figure  D. 

b.  The  Least  Squares  Regression  Procedure  to  Solve  for  {#,.} 

The  previous  estimates  of  rf}  {pj}  and  {$?}  are  input  to  this  stage. 
The  partial  derivative  of  equation  (3.2)  with  respect  to  4  is 


<31nL 

ap* 


ainir,.  c(.+^°  as/ 

S^e"0 


(3.16) 


where 
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Thus 


M  Ni 


fr££'l 

dP*  b  1  >1.  I  5i+c.7>  LV  v  7 


Ivi  •  (3-i7> 


Two  terms  of  a  Taylor  expansion  of  equation  (3.17)  yields 


a,  T  m  Ni 

o-^-Z  E 


ap 


k  M  J-l 


-(exp] 


(yr2/)e'??fe  v^'K-P?  V 

v  JA*0 
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)• 


(3.18) 


Let 


= 


1 


(3.19) 


The  equation  (3.18)  can  be  rewritten  as 
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21  T  M  N> 

0=SEt*YY 
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h*  0 


(3.20) 


Let 


uiih  =  xah»nc 


< 


(3.21) 


and 


(3.22) 


The  equation  (3.20)  can  be  rewritten  as 


M  Ni  (  p 

0  =  EE  za  -  Em(,aPa 


M  /« 1 


A*  0 


% 


(3.23) 


which  are  the  normal  equations  for  a  least  squares  regression  having  dependent 
variable  and  independent  variables  Up,  k=0,l>->p. 

The  solution  of  (3.23)  is 


B  =  [TJJ)'XITZ 


(3.24) 


where  Z  and  U  are  defined  as  in  equations  (3.25)  and  (3.26)  respectively. 
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This  procedure  is  an  iteratively  re-weighted  least  squares  regression;  cf.  McCullagh 
and  Nelder  [Ref.  7]. 


c.  Newton  Procedure  to  Solve  for  {£;) 

The  prior  estimates  of  rf,  and  are  input  to  this 

stage.  In  order  to  update  the  values  of  {^},  we  need  to  have  the  partial  derivative 


of  unconditional  ln-likelihood  function  for  observer  i  with  respect  to  The 


unconditional  ln-likelihood  function  for  observer  i  is  the  same  as  equation  (2.11). 
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The  derivative  is  as  follows : 


aih i. 


4 


35.. 


U=i 

n> 


C,+eT)  35, 


5,+e’1° 


i  y-i 


(3.27) 


Again,  let 


ou-Si 


rv  *  (p*  -  */)e 


The  equation  (3.27)  can  be  rewritten  as 


w, 


31nl.  _ 

-aT  -  »«,>  -  -c«  +  E>V 

;»1 


C.+e^ 

-V7— ;exP(r«> 

SAe' 


(3.28) 


—  „  u 

C,+e^ 

Treating  — -  as  a  constant,  then 


S.+e" 


&kxL.  dh(U  2 

—rm*-Ci  -  -1— (V1? 

3$?  5.+C*'  /*» 


(3.29) 


31ni. 

If  - 1  =  h(Q  =  0  ,  the  Newton  equation  for  ^  is 

3^. 


31nJL. 


0.= -*(?,)  = 


(3.30) 
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Update  the  value  of  £■  by  solving  the  following  equation  : 


(3.31) 


As  noted  before,  it  is  possible  for  equation  (3.31)  to  be  badly  conditioned.  If 

A 

—hkfj  <  e2  ,  where  1  >  >  c1>e2>0;  then  equation  (3.31)  is  badly 

behavior  happens  occasionally  and  occurs  for  only  one  £j  out  of 

d  ,/,.o\| 


?)<*! 


and 

conditioned.  This 
the  M.  To  avoid  this  behavior,  if 


?)<ei 


and 


<e2  ,  then  put 


£,=  5,  ;  otherwise  let  Ci.be  the  solution  to  equation  (3.31).  Return  to  the  bisection 
search  procedure  and  iterate  steps  3.a  to  3.c  until 


/ 

max 

\ 


i 

> 

o  — 
yy* 

/ 

<  0.01  for  all  k=0,..,p  and  i  =  1,..,M. 
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C.  IMPLEMENTATION  OF  SIMULATION 

Simulation  will  be  used  to  study  the  small  sample  behavior  of  the  estimators 
of  the. five  parameters  of  the  gamma/Weibull  model.  Different  numbers  of  observers 
(M)j  targets  (NJ  and  opportunity  times  (0^)  are  used  as  the  input  to  the  simulations. 
For  each  set  of  input,  (0,M,Nj),  the  simulation  run  starts  with  the  same  random  seed 
as  in  Chapter  II.  The  simulation  is  run  until  there  are  100  (R)  estimates  of  r?.  A  run 
will  contain  more  than  100  replications  if  the  simple  Weibull  regression  model  is 
found  to  describe  the  variability  of  the  data  for  some  of  the  replications  or  if  the 
procedure  to  find  the  estimate  of  r\  is  badly  conditioned.  The  statistics  of  mean  bias 
(M.B),  mean  square  error  (M.S.E)  and  their  respective  standard  errors  (  S.E(M.B) 
and  S.E(M.S.E)  )  for  all  the  estimates  are  computed  using  the  definitions  in 
equations  (2.16a,b)  and  (2.17a,b).  The  averaged  uncensoring  levels  (UC)  as 
computed  using  equation  (2.18)  are  also  reported.  The  statistics  for  {Q  do  not 
include  those  ^  for  which  |h(^)|  <  e„  (e, =0.1),  and  |h'(£j)|  <  (e2=0.01). 

The  remainder  of  this  section  provides  an  outline  of  one  replication  in  the 
simulation. 

1.  Input  to  Simulation 


•  M  the  number  of  observers  (  Possible  values  are  15 , 25  and  35  ). 

»  £j  the  shape  parameter  of  the  Weibull  distribution  for  ith  observer, 

where  i  =  1,2,. ..M  ;  (  Set  =  0  for  all  i ). 

•  rj  the  parameter  for  the  single-parameter  gamma  distribution, 

GAM(a,o) ,  where  a  =  en  ;  (Set  r?  =  1 ). 

•  N;  the  number  of  targets  presented  to  the  ith  observer  ,  where  i  = 

1,2,...,M  ;  (  Possible  values  are  15  ,25  and  35). 
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•  xijk  the  kth  environmental  factor  which  affects  the  observation  time 

of  the  ith  observer  when  presented  with  the  jth  target. 

•  Ojj  the  opportunity  time  for  the  ith  observer  to  detect  the  jlh  target; 

(Possible  values  are  10  and  15). 


•  {/}k}  regression  coefficients  for  the  scale  parameter  of  Weibull 
distribution  ,WEI(  ,  e'li  ),  where  k  =  0,l,2,...,p  and 


p..  =  expEx^pJ 

\k*  0  J 

=  0.5  ). 


;  ( set  p  =  2  ;  =  0.8  ,  ^  =  -0.2  and  p2 


Simulation  of  Data  for  the  i‘h  Observer 


•  Generate  a  single-parameter  gamma  r.v.  Aj  from  GAM(a, a), having 
density  function  as  in  equation  (2.1)  where  a  =  e11  and  rj  =  1. 

•  Generate  environmental  factors  r.v.  xijk ;  the  Xyk  are  independent  with 
normal  distributions  having  mean  /ik  and  variance  o\  for  k  = 

1,2,. ..,p;  (Set  p  =  2  ;  Pj  =  1  ,  Oj  =0.5  and  =  2  ,  o\  =  1). 
Put  xij0  =  1  for  all  i,j. 

•  Compute  the  Weibull  shape  parameter  as  in  equation  (3.1). 

•  Generate  independent  exponential  distribution  r.v.’s  with  mean 
1;  i.e.  W;j  ~  EXP(l)  for  all  i  and  j. 


Compute  the  target  detection  times :  1C  = 


\Ai  / 


•  Compute  the  recorded  ln-times  Yy  ,  censoring  indicator  Ay  and  Q 
using  the  definitions  in  equation  (2.3) ,  (2.4)  and  (2.5)  respectively. 
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3.  Initial  Estimate  of  £k’s  and  £{  in  the  Simple  Weibull  Regression  Model 
without  the  Hierarchical  Gamma  r.v.’s 


Step-1.1  :  Initially,  set  =  In U{j  and  =0 


Step-1.2  :  Iteration  for  regression  : 


•  Compute  tjj ,  Wjj ,  uSj  and  zti  as  in  equations  (3.3)  to  (3.6). 

•  Regress  the  dependent  variable  zYl  against  the  independent 
variables  uijh  for  {/JJ  as  in  equations  (3.7)  to  (3.9). 

Step-1.3  :  Newton’s  procedure  for  estimating 

•  Update  the  value  of  r^  as  in  equation  (3.10). 

•  Compute  the  new  as  in  equation  (3.11)  with  no  need  to  check 
the  ill  condition. 

Step-1.4  Put  =  {!-.}  and  {p°}  =  for  all  i  =  l,2...,M  and 

p 

k = 0,1, ..,p.  Compute  new  value  of  2„£°  =  T,  xijk  p*  and  return 
to  Step- 1.2  for  only  one  iteration. 

Step-1.5  Check  for  extra-variability  in  the  data  by  computing  the  sample 
variance  d\  for  the  times  Uy  as  in  equation  (3.12). 

If  bE  <  1  ,  stop  and  quote  the  simple  Weibull  regression 
model  only  with  the  parameters  as  estimated. 
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Otherwise,,  execute  next  the  recursive  procedure. 

4.  Recursive  Procedure  to  Estimate  the  Parameters  of  the  Full  Hierarchical 
Gamma/Weibull  Model 


Input :  initial  p*  ,k=0,l,...,p  ;  ,i=l,2,....,M. 

Step-2.1  :  Initial  condition  for  gamma  parameter  r\ 

•  Update  x  B°  =  £x„fcpj  . 

v  *■  o 

•  Find  the  initial  condition  for  r?  using  equations  (3.13)  and  (3.14) 

If  <0  ,  set  r?°  =  0. 

6 

Step-2.2  Perform  the  bisection  search  for  rj  using  the  algorithm  shown  in 
Figure  D  (next  page).  The  algorithm  includes  a  check  for  the 
flatness  of  the  curve  of  f (rj). 
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Figure  D  The  Algorithm  of  Bisection  Search  with  Instability  Check  for  t) 


11  In  the  Algorithm,  \vc  let  // 


1)) 


and  k  is  the  current  index  of  increments. 


6 


,  =n0+*x$ 

where  u  0  „  . 

tiw-irT +(*-1)x6 


Step-3.1  :  Iteration  for  the  regression  coefficients  ,{/3k) 

•  Compute  w-j  ,  Ujj  and  zVl  using  equations  (3.19),  (3.21)  and 
(3.22). 

•  Regress  the  dependent  variable  Zjj  against  the  independent 
variables  uijh  using  equations  (3.24) ,  (3.25)  and  (3.26). 

Step-4.1  ;  Iteration  for  the  ln-shape  parameters  ,{£;},  of  the  Weibull' 
distribution 


Update  xfl  =  Zxtfl  and  ri}  =  fa  - 


•g 

0 


Compute  )  and 
(3.29)  respectively. 


_9  dr® 

di 


using  equations  (3.28)  and 


•  If  /»(£,•)<  Ej  and 


<e2  ,  for  =  0.1  and  e,  -- 


0.01  ,  keep  .  Otherwise,  update  using  equation 

(3.31). 

Step-4.2  :  Test  for  stop  criteria 


/ 


If  max 


ni-i h 


9a/ 


P*-Pi 


P! 


<  0.01  for  all  k=0,..,p  and 


i  =  l,..,M,  then,  stop  and  record  the  estimates.  Otherwise,  set 
T)°  =  ri  »  Pt  =  Pjt  >  and  return  to  Step-2.2. 
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The  simulation  programs  for  this  gamma/Weibull  regression  model,  named’ 
"SIMULA2",  "GAMWEI"  and  "FVALUE",  are  detailed  in  Appendix  D. 

D.  RESULTS 

Results  of  the  simulation  experiments  are  presented  in  Appendix  B.  The  results 
are  presented  in  five  sections.  Each  section  shows  results  for  one  estimator;  the 
gamma  parameter  i)  and  the  Weibull  parameters  0O ,  fa  ,  ft2  and  £;.  Each  section 
contains  tables  of  statistics  and  histograms  of  estimate  bias  for  different  numbers  of 
observers  (M  or  OBS),  targets  (N  or  TGT)  and  opportunity  times  (O). 

The  first  two  tables  in  each  section  present  the  statistics  for  different 
opportunity  times  ( 0= 10  and  0= 15  in  our  case ).  The  next  two  tables  organize  the 
information  to  highlight  relations  between  the  number  of  targets  and  the  number  of 
observers  and  the  statistics  of  the  estimates.  Histograms  of  the  biases' of  the 
estimates  are  also  presented.  All  parameter  estimates  appearing  in  Appendix  B  are 
from  replications  for  which  the  full  hierarchical  gamma/Weibull  is  estimated. 
Since  5,=1  ,  the  statistics  for  {^}  are  for  the  average  of  the  {fj  for  all  the 
observers  in  a  replication.  The  average  includes  only  those  £,’s  for  which  the  Newton 
procedure  converged.  The  replications  for  which  62E  <  1  and  the  replications  in 
which  rj  could  not  be  found  are  not  used  to  compute  the  statistics  appearing  in 
Appendix  B;  statistics  concerning  the  numbers  of  these  cases  appear  in  Tables  A  and 
B  below. 
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The.statisticsin  Appendix  B  use  100  replications.  Tables  A  and  B  below  show 
-  the  averaged  uncensoring  levels,,  the  number  of  additional  replications  for  which  rj 
could,  not  be  found  due  to  the  flatness. of  f,  and  the  number  of  times  a  £j  could  not 
be  found  using  Newton  procedure  due  to  the  smallness  of  |  h(^)  |  and  |  h'  (£;)  | .  The 
results  of  Tables  A  and  B  (next  page)  are  summarized  as  follows : 


•  Uncensoring  level  (UC) : 

A  longer  opportunity  time  (O)  is  reflected  in  a  higher  average  uncensoring 
level  (UC).  The  opportunity  time  O  =  10  gives  about  76%  UC,  while  O 
=  15  gives  about  85%  UC.  This  behavior  is  what  we  would  expect. 

•  The  variability  of  the  data  explained  by  the  simple  Weibull  regression  : 

Different  parameter  estimates  in  the  simple  Weibull  regression  cause  the 
different  number  of  cases  for  which  {  b\  <  1}  for  the  same  input  of 
number  of  observers  and  number  of  targets  in  Tables  A  and  B.  In  both 
Tables  A  and  B,  the  cases  of  {  b\  <1}  appear  more  frequently  when 
the  number  of  observers  (M)  is  smaller.  This  is  reasonable  behavior  since 
more  observers  tend  to  provide  more  evidence  that  a  hierarchical  model 
is  needed  to  explain  the  variability  of  the  data.  The  manner  in  which 
6 ^  is  computed  prevents  drawing  conclusions  concerning  the  effect  of 
censoring  on  the  assessment  of  unexplained  variability. 

•  Instability  check  for  r? : 

There  are  only  two  replications  in  which  rj  could  not  be  found  using  the 
bisection  search  for  the  root  of  f(r?).  The  function  f(rj)  for  one  of  them 
is  plotted  as  in  Figure  C,  and  the  function  for  the  other  is  similar.  Both 
replications  used  the  same  random  seed  but  had  a  different  set  of  inputs, 
(M,Nj,Ojj).  These  two  replications  were  stopped  and  new  random  numbers 
were  drawn.  This  behavior  did  not  occur  in  the  model  of  Chapter  II  with 
known  parameters  for  the  Weibull  distribution. 

•  Instability  check  for  £  : 

Smaller  opportunity  times  tend  to  have  more  replications  in  which  the 
procedure  to  find  £j  is  badly  conditioned.  The  maximum  number  for  any 
case  is  2. 
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TABLE  A :  Comparisons  of  Uncensoring;  levels,  Extra  Variability  Checks  and  Checks 


Number  of 
,  Observers 

Number  of 
Targets 

Average 

ud 

Number  of 

<  *i  <  i> 

Num.  unstable 

m 

Num.  unstable 
h  and  h'  12 

15 

.15 

76;4% 

li 

0 

.  la 

■Em 

m 

5  . 

1 

• .....  0 

35 

75.6% 

6 

0 

1 

25 

15 

75.8% 

1 

0 

lb 

25 

75.8% 

2 

0 

1c 

35 

76.4% 

1 

0 

0 

35 

15 

76.2% 

0 

0 

2d 

.  25 

76.3% 

1 

le 

35 

76.3% 

0 

0 

0 

TABLE  B  :  Comparisons  oiTJncensoring  levels,  Extra  Variability  Checks  and.Checks- 


Number  of 
Observers 

Number  of 
Targets 

Average 

UC 

Number  of 

{ <  i> 

Num.  unstable 

m 

Num.  unstable 
h  and  h' 

15 

15 

85.0% 

23 

0 

0 

25 

85.4% 

6 

1 

0 

35 

84.6% 

11 

0 

0 

25 

15 

84.5% 

3 

0 

If 

25 

84.7% 

7 

0 

0 

35 

85.3% 

5 

0 

0 

35 

15 

85.2% 

5 

0 

0 

25 

84.9% 

3 

0 

0 

35 

85.1% 

2 

0 

0 

1 

"  All  the  cases  of  numerically  unstable  h  and  h'  were  obtained  while  replicating  the  simulations 
starting  with  random  seeds  a:999445582,  b:444977940,  c:5387909S6,  d:999445582  &.  153499960, 
e:1011786849  and  f:444977940. 
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In  Appendix  B,  the  statistics  for  the  average  observer  ^  do  not  include  those 
individual  £;  for  which  the  Newton  procedure  did  not  converge.  The  statistics  for  pk 
do  not  include  those  replications  for  which  the  simple  Weibull  regression  model  was 
used  to  describe  the  data.  A  summary  of  the  results  of  the  simulation  experiment 
shown  in  Appendix  B  follows : 


•  The  estimate  of  the  gamma  parameter,  fj 

The  mean  bias  tends  to  be  positive.  The  histograms  of  biases  indicate  that 
the  distribution  is  skewed  to  the  right.  Since  the  variance  of  A  is  e‘\  the 
positive  mean  bias  of  r?  suggests  that  for  small  sample  sizes  the  procedure 
is  indicating  less  between  variability  than  there  is  in  the  data.  Increasing 
the  number  of  observers  tends  to  decrease  the  mean  bias  and  mean  square 
error.  Increasing  the  opportunity  time  has  little  effect  on  mean  bias  and 
mean  square  error.  Increasing  the  number  of  targets  has  little  effect  on 
mean  bias  and  mean  square  error.  As  noted  in  Chapter  II,  this  behavior 
is  reasonable. 


•  The  estimates  of  the  Weibull  parameters,  p0  through  p2 

The  first  two  tables  in  parts  2  through  4  of  Appendix  B  indicate  that  all 
the  M.B  and  M.S.E  of  {/1J  are  quite  small  compared  to  those  for  rj.  This 
is  due  to  the  larger  sample  size  used  in  computing  {jSJ.  An  increase  in 
the  opportumty  time  has  little  effect  on  the  mean  bias  and  mean  square 
error.  Increasing  the  number  of  targets  or  observers  with  a  fixed  number 
of  the  other  yields  small  changes  in  M.B  and  M.S.E.  However,  increasing 
both  the  number  of  targets  and  the  number  of  observers  does  tend  to 
decrease  the  mean  bias  and  mean  square ‘error  of  {/JJ.  The  M.B’s  are 
both  positive  and  negative. 


1  M  * 

•  The  average  estimate  of  the  Weibull  parameter,  £=— 

Mm 

Since  =1  in  the  simulations,  the  histograms  and  statistics  include  the 
averaged  for  all  the  observers  for  which  the  Newton  procedure  for  fj 
converged.  Generally,  all  the  magnitudes  of  M.B  and  M.S.E  of  l  in  the 
first  2  tables  of  part  5  in  Appendix  B  are  quite  small.  All  the  M.B  of  the 
averaged  1  are  negative.  The  M.B.  and  M.S.E.  are  slightly  smaller  for 
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longer  opportunity  times.  Increasing  the  number  of  targets  or  the  number 
of  observers  with  the  other  fixed  results  small  changes  in  the  M.B  and 
M.S.E,  Increasing  both  the  number  of  targets  and  the  number  of 
observers  tends  to  decrease  the  mean  bias  and  mean  square  error. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  main  effort  of  this  thesis  is  to  study  the  small  sample  behavior  of 
estimators  for  a  hierarchical  gamma/Weibull  regression  model  for  target  detection 
times,  this  model  can  also  be  used  to  describe  and  predict  times  to  failure  of  similar 
machines,  (e.g.  engines)  in  different  environments.  In  this  thesis,  the  model  assumes 
that  there  are  M  observers.  The  i,h  observer  is  presented  with  Nj  targets.  The 
variable  Uy  is  the  time  it  takes  for  observer  i  to  detect  the  target  j.  Let  xijl ,  xij2 ,...., 
Xjjp  be  the  values  of  environmental  variables  which  may  affect  Uir  Given  the 
quantified  value  of  the  ability  for  observer  i,  Aj  =  0,  the  Ujj’s  are  assumed 
conditionally  independent  random  variables  having  Weibull  distributions, 
WE^,e~^\  ,  where  the  scale  parameter  has  the  form 
=  expji^jij  =  expjEx^p*!  .  The  random  variables  {A}  are  assumed 
independent  from  a  single  parameter  gamma  distribution ,  GAM(a,a)  with  a  -  e\  A 
numerical  procedure  based  on  maximum  likelihood  is  used  to  estimate  the 
parameters  of  the  model.  The  numerical  procedure  is  iterative  and  uses  a  bisection 
root  search  method  for  estimating  the  gamma  parameter  r? ,  least  squares  regression 
to  estimate  the  Weibull  In-scale  parameters  P0  through  p2  and  an  approximate 
Newton  procedure  to  estimate  the  Weibull  ln-shape  parameters  {£,}.  Simulation  is 
used  to  study  the  behavior  of  the  estimates  for  small  sample  sizes.  Generally,  the 
numerical  estimation  procedures  work  well.  There  are  few  numerical  problems.  All 
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the  mean  biases  and  mean  square  errors  of  the  estimates  tend  to  be  small  compared 
to  the  theoretical  inputs  to  the  simulation.. 

Topics  for  future  research  include  replacing  the  one  parameter  gamma 
distribution  with  the  two  parameter  gamma  distribution  in  the  hierarchical  model  and 
exploring  the  use  of  EM  methodology  to  estimate  model  parameters.  The  two 
parameter  gamma  distribution  arises  as  the  posterior  distribution  of  the  current 
model.  Another  area  is  to  study  using  the  fitted  hierarchical  gamma/Weibull  model 
to  predict  future  performance. 

It  is  hoped  that  the  hierarchical  gamma/Weibull  regression  model  will  be  a 
useful  tool  to  describe  and  predict  one  aspect  of  the  effect  of  human  performance 
on  the  battlefield. 


43 


LIST  OF  REFERENCES 


1.  Gaver,  D.  P.,  Jacobs,  P.  A.,  and  O’Muircheartaigh,  I.,  Hierarchical  Models  for 
Target  Detection  Times,  Naval  Postgraduate  School  Technical  Report.  To 
appear. 

2.  Bae,  Deok  Hwan,  Models  for  Target  Detection  Times,  Naval  Postgraduate  School 
Thesis,  September  1989. 

3.  Johnson, L.  W.  and  -Riess,  R.  D.,  Numerical  Analysis,  Addison-Wesley,Inc.  New 
York,  1977. 

4.  Peter  A.  Stark, Introduction  to  Numerical  methods,  Macmillan  Publishing  Co., Inc. 
New  York  ,1970. 

5.  Lewis,  P.  A.  and  Uribe, L.,77ie  New  Naval  Postgraduate  School  Random  Number 
Package  LLRANDOM II,  Naval  Postgraduate  School  Technical  Report,  NPS55- 
81-005,  Monterey,  CA,  1981. 

6.  Leonard  Gilman  and  Allen  J.  Rose,  APLAn  Interactive Approach,Third Edition, 
JOHN  WILEY  &  SONS, Inc.  New  York,  1984. 

7.  P.  McCullagh  and  J.A.  Nelder,  Generalized  Linear  Models,  Chapman  and  Hall, 
New  York,  1983. 


44 


APPENDIX  A.  STATISTICS  AND  GRAPHICS  FOR  THE  ESTIMATE  OF 
>/:  IN  THE  SINGLE  PARAMETER  GAMMA  MODEL 


Table  1.  Mean  Bias,  Mean  Square  Error,  Std.  Error  and  Uncensored  Level  for  tj  at 
6  =  10  in  the  Gamma  Model 


Number  of 

■SI 

A 

'/ 

uc 

Observers 

M.S.E(S.E) 

15 

0.0727  (0.0497) 

0.2477  (0.0693) 

81% 

15 

25  , 

0.1027  (0.0436) 

0.1967  (0.0269) 

82% 

35 

0.067  (0.0466) 

0.2173  (0.0569) 

82% 

.15 

0.1003  (0.0311) 

0.1051  (0.0167) 

S2% 

25 

25 

0.0566  (0.0304) 

0.0935  (0.012  ) 

82% 

35 

0.06S2  (0.0274) 

0.0784  (0.0108) 

82% 

15 

0.058  (0.0281) 

0.0S07  (0.0129) 

82% 

35 

HM 

0.00SS  (0.0253) 

0.0626  (0.0092) 

S2% 

35 

0.0472  (0.0221) 

0.05  (0.0103) 

82% 

Table  2.  Mean  Bias,  Mean  Square  Error,  Std.  Error  and  Uncensored  Level  for  ;/  at 
0=15  in  the  Gamma  Model 


Number  of 
Observers 

Number  of 
Targets 

A 

'/ 

UC 

M.B(S.E) 

M.S.E(S.E) 

15 

15 

0.0614  (0.0434) 

0.1886  (0.0413) 

89% 

25 

0.0884  (0.0428) 

0.1875  (0.026  ) 

90% 

35 

0.053  (0.0452) 

0.2032  (0.052  ) 

90% 

25 

15 

0.094  (0.0305) 

0.0998  (0.0151) 

90% 

25 

0.0517  (0.03  ) 

0.0911  (0.0115) 

90% 

35 

0.0706  (0.0274) 

0.0785  (0.0103) 

90% 

35 

15 

0.0519  (0.0277) 

0.0781  (0.0132) 

90% 

25 

0.0085  (0.0246) 

0.0593  (0.0086) 

90’% 

35 

0.0471  (0.0216) 

0.0481  (0.009  ) 

90% 

Table  3.  Tendencies  of  Mean  Bias  and  Mean  Square  Error  for  //  at  O  =  1G  in  the 
Gamma  Model, 


M.B(M.S.E)  of  // 

15  TGT 

25  TGT 

35  TGT 

'  15  OBS 

'  0.d73(Q.24S) 

0.103(0.197) 

0.067(0.217) 

25  OBS 

0:1,(0.105) 

0.057(0.094) 

0.068(0.078) 

35  OBS 

0.058(0.081) 

0.009(0.063) 

0.047(0.05  ) 

Figure  1.  The  Tendencies  of  the  Bias  of  ?/  with  changes  in  OBS  and  TGT  at  0=  10 
in  the  Gamma  Model:  M  =  number  of  observers  ;  X  =  number  of 
targets. 
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N  =15  N  =  25  N  =  35  :  M  =  15  M  =  25  M  =  35 

Figure  2;  The  Tendencies  of  the  Bias  of  ?/  with  changes  in  OBS  and  TGT  at  0=  10 
in  the  Gamma  Model:  M  =  number  of  observers  ;  K  =  number  of 
targets. 


Figure  3.  Tiie  Tendencies  of  the  Bias  of ;/  with  changes  in  OBS  and  TGT  at  O  =  10 
in  the  Gamma  Mode!:  M  =  number  of  observers  ;  N  =  number  of 
targets. 


TGI  =  15  ,  0  =  10  (  UC  :  82  PERCENT  ) 


Figure  4.  Comparison  of  the  Bias  of  ?/  between  Different  numbers  of  Observers  with 
TGT  =  15  at  0=10  in  the  Gamma  Model:  TGT  =  the  number  of 
targets. 
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Figure  5.  Comparison  of  the  Bias  of  >/  between  Different  numbers  of  Observers  with 
TGT  =  25  at  0=10  in  the  Gamma  Model:  TGT  =  the  number  of 
targeis. 


V' 


Figure  6.  Comparison  of  the  Bias  of ;/  between  Different  numbers  of  Observers  with 
TGT  =  35  at  0=  10  in  the  Gamma  Model:  TGT  =  the  number  of 
targets. 
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Table  4.  Tendencies,  of  Mean  Bias  and  Mean  Square  Error  for  >/  at  O  =  15  in  the 
Gamma  Model _ _ '  _ _ 


j  M.B(M.S.E)  of )/ 

15  TGT 

25  TGT 

35  TGT 

15  OBS 

0.061(0.189), 

0.088(0.188) 

0.053(0.203) 

25 OBS 

0.094(0.1  ) 

0.052(0.091) 

0.071(0.078) 

35  OBS 

0.052(0.078) . 

0.008(0.059) 

0.047(0.04$) 

0  =  15  (  UC  :  60  PERCENT  ) 

M  -  15  <=>  N  =  15 


* 


N  =  15  N  =  25  N  o  35  :  M  »  15  M  =  25  M  =  35 


Figure  7.  The  Tendencies  of  the  Bias  of ;/  with  changes  in  OBS  and  TGT  at  0=15 
in  the  Gamma  Model:  M  =  number  of  observers  ;  X  =  number  of 
targets. 
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C!\S  OT  ETA  fi  I  G'AS 


Figure  10.  Comparison  of  the  Bias  of  ?/  between  Different  numbers  of  Observers 
with  TGT  =  15  at  0=  15  in  the  Gamma  Model:  TGT  =  the  number 
of  targets, 
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Figure  II.  Comparison  of  the  Bias  of  ij  between  Different  numbers  of  Observers 
With  TGT  =  25  at  0=  15  in  the  Gamma  Model:  TGT  =  the  number 
of  targets. 
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APPENDIX  B.  STATISTICS  AND  GRAPHICS  FOR  ALL  THE 
ESTIMATES  IN  GAMMA/ WEI  BULL  REGRESSION  MODEL 

1.  TABLES  AND  HISTOGRAMS  FOR  >/ 


Table  5.  Mean  Bias,  Mean  Square  Error  and  Std.  Error  for  tj  at  O  =  10  in  the 
_ GAM /WEI  Regression  Model  _ ' _ _ 


Number  of. 

(  Observers 

Number  of 
Targets 

A 

■V 

M.B(S.E) 

M.S.E(S.E) 

H 

15 

.1770(.0552)  , 

.3299(!0S92) 

25 

,1573(.0455) 

.2274(.0445), 

35 

.0761  (.0440) 

.  1 955(.0290) 

25 

15 

.0947(.0419) 

.  1 S 1 1  (.0420) 

25 

,0571(,0315) 

.1003(.0206) 

35 

,1422(.0324) 

.1228(.0208) 

35 

15 

.1074(.0352) 

.  1 32S(.02 1 6) 

25 

.0658(.0307) 

,0968(,0129) 

35 

.0940(.0223) 

.0577(.0070) 

Table  (>,  Mean  Bias,  Mean  Square  Error  and  Std.  Error  for  ;/  at  0  =  15  in  the 
_  GAM/ WEI  Regression  Model  _ _ 


Number  of 
Observers 

Number  of 
Targets 

A 

>/ 

M.B(S.E) 

M.S.E(S.E) 

15 

15 

,2107(.0561) 

.3524(.1023) 

25 

.1934(.0451) 

.2365(.045S) 

35 

.11 17(.0442) 

.2038(.0319) 

25 

15 

.1 170(.0418) 

.1851  (.0454) 

25 

.0755(.0314) 

.  102 1(.0222 ) 

35 

.1434(.0314) 

.1 1 74(.0 1 97) 

35 

15 

.1450(.0336) 

.1315(.0228) 

25 

.0800(.0306) 

.0982(.0133) 

35 

.0968(.0219) 

.0566(.0072) 
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Figure  13.  Comparison  of  the  Bias  of ;/  between  Different  numbers  of  Observers 
with  TGT=  15  at  0=  10  in  the  GAM/WEI  Regression  Model:  TGT 
=  the  number  of  targets. 


Figure  14.  Comparison  of  the  Bias  of  >j  between  Different  numbers  of  Observers 
with  TGT=  25  at  0=  10  in  the  GAM/WEI  Regression  Model:  TGT 
=  the  number  of  targets. 
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Figure  15.  Comparison  of  the  Bias  of  >/  between  Different  numbers  of  Observers 
with  TGT=35  at  0=  10  in  the  GAM/WEI  Regression  Model:  TGT 
=  the  number  of  targets. 


Table  8.  Tendencies  of  Mean  Bias  and  Mean:  Square  Error  for  ?/  at  O  =  15  in  the 


QAM /WEI  Regression  Model 


Figure  16.  Comparison  of  the  Bias  of  >/  between  Different  numbers  of  Obseners 
with  TGT=  15  at  0=  15  in  the  GAM/WEI  Regression  Model:  TGT 
=  the  number  of  targets. 
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Figure  17.  Comparison  of  the  Bias  of  >/  between  Different  numbers  of  Obseners 
with  TGT=25  at  0=  15  in  the  GAM/WEI  Regression  Model:  TGT 
=  the  number  of  targets. 
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Figure  18.  Comparison  of  the  Bias  of »/  between  Different  numbers  of  Observers 
with  TGT=35  at  0=  15  in  the  GAM/WEI  Regression  Mode!:  TGT 
=  the  number  of  targets. 
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2.  TABLES  AND  HISTOGRAMS  FOR  /JO 


Table  9.  Mean  Bias, Mean  Square  Error  and  Std.  Error  for  ,/10  at  0  =  10  in  the 
GAM/ WEI  Regression  Model 


•  Number  of 

Number  of 

A 

k 

Observers 

Targets 

M.B(S.E) 

M.S.E(S.E) 

15 

-.0605(.0317)  ‘ 

,1021(.0123) 

15 

25 

„  -.0234(.0278) 

.076400109) 

35 

.0304(.0223> 

.049900067) 

15 

i0458(.023S) 

.0575000S1) 

25 

25 

-;0167(.0222) 

,0487(.0052) 

35 

-.009400197) 

.038300056) 

15 

-.0067(.019S) 

.03S6(.0067) 

;  35 

25 

-.039800176) 

.0319(.0047) 

35 

.0062{.0152) 

.022700027) 

Table  10.  Mean  Bias, Mean  Square  Error  and  Std.  Error  for  /JO  at  O  =  15  in  the 
,  GAM/WEI  Regression  Model 


Number  of 

Number  of 

A 

k 

Observers 

Targets 

M.B(S.E) 

M.S.E(S.E) 

15 

-.060900316) 

.101500127) 

15 

25 

-.023400271) 

.072500095) 

35 

.026300224) 

.049800063) 

15 

.0365(.0244) 

.059500073) 

25 

25 

-.029200210) 

.043900049) 

35 

-.01S5001S4) 

.03351.0053) 

-.004300191) 

.03581.0066) 

35 

25 

-.045400176) 

.03241.0041) 

35 

-.002100138) 

.018700023) 

Table  11.  Tendencies  of.Mean  Bias  ahd;  Mean  Square  Error  for /?0  at  O  =  10  in  the 
GAM/WEI  Regression  Model 


M.B(M.S.E)  of  ji0.  , 

15  TGT 

25  TGT 

35  TGT 

1;5  OBS 

-.061  (.102) 

-.023(.076) 

.030(.050) 

25  OBS 

.046(.05/V 

-.017(.049), 

-.009(.038) 

35  OBS 

-,007(.039) 

-;040(.032) 

.006(.023) 

Figure  19.  Comparison  of  the  Bias  of  p0  between  Different  numbers  of  targets  with 
OBS=  15  at  0=  10  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 
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Figure  20.  Comparison  of  the  Bias  of  /JO, between  Different  numbers  of  targets  with 
OBS  =  25  at  0=  10  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 


Figure  21.  Comparison  of  the  Bias  of  fiO  between  Different  numbers  of  targets  with 
OBS  =  35  at  0=  10  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 
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Table  12.  Tendencies  of  Mean  Bias  and  Mean  Square  Error  for  /JO  at  O  =  15  in  the 
_  GAM/ WEI  Regression  Model  _ _ 


M.B(M.S.E)  of  p0 

15  TGT 

25  TGT 

35  TGT 

15  OBS 

-.061 (.102) 

-.023(.073) 

,026(.050) 

25  OBS- 

.0370060) 

-.029(.044) 

-.01 8(.033) 

35  OBS 

-.004(.036) 

-.045(.032) 

-.002(,019) 

«b 


z 


OBS  =  15 , 0  = 
N  -  15 


15  (  UC  :  85  PERCENT  ) 

N  -  25 


Figure  22.  Comparison  of  the  Bias  of  /JO  between  Different  numbers  of  targets  with 
OBS=  15  at  0=  15  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 
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Figure  24.  Comparison  of  the  Bias  of  fiO  between  Different  numbers  of  targets  with 
OBS  =  35  at  0=  15  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 


3.  TABLES. AND  HISTOGRAMS  FOR  /ll 


Table  13.  Mean  Bias, Mean  Square  Error  and  Std.  Frror  for  pi  at  O  =  10  in  the 
_  GAM/WEI  Regression  Model 


Number  of 

Number  of 

A 

Pi 

Observers 

Targets 

M.B(S.E) 

M.S.E(S.E) 

15 

.0180(.0157) 

.0243(.0032) 

15 

25 

-.0002(.0140) 

.0193(.0026) 

35 

-.0010(.0112) 

.0122(.0017) 

15 

-,0154(.0132) 

,0173(.0026) 

25 

25 

-.0006(.0096) 

.0091 (.00 14) 

35 

-0036(.0080) 

.0063(.0010) 

15 

.0188(.0126) 

.0160(.0020) 

35 

25 

.0094(.0096) 

,0091  (.0011) 

35 

.0006(.0069) 

.0047(.0006) 

Table  14.  Mean  Bias,Mean  Square  Error  and  Std.  Error  for  pi  at  O  =  15  in  the 
_ GAM/WEI  Regression  Model _ 


Number  of 

Number  of 

A 

Pi 

Observers 

Targets 

M.B(S.E) 

M.S.E(S.E) 

15 

.0079(.0145) 

.0205(.002S) 

15 

25 

-.0060(.0127) 

.0158(.0022) 

35 

,0076(.0107) 

.01 13(.0014) 

15 

-.0214(.0127) 

.0162(.0024) 

25 

25 

.0036(.0090) 

.0079(.0013) 

35 

-.0046(.0075) 

.0055(.0009) 

15 

.0150(.0121) 

.0146(.0018) 

35 

25 

.0077(.0092) 

.0084(.0010) 

35 

.0009(.0067) 

.0043(.0005) 
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Table  15.  Tendencies  of  Mean  Bias  and  Mean  Square  Error  for  fll  at  O 
GAM/WEI  Regression  Model 


=  10  in  the 


15  OBS 
25  OBS 
35  OBS 


15  TGT 

,bl8(.024) 

-.015(.017) 

.019(1016) 


M.B(M.S.E)  of  /?, 
25  TGT 
■0Qb(.Q19) 
•■001(.009.) 
.009(.009) 


35  TGT 
-.00K.012) 
-,004(.006) 
.001(.005) 


OBS  = 
N  -  15 


C 

BIAS  or  t 


Figure  25.  Comparison  of  the  Bias  of  0 1  between  Different  numbers  of  targets  with 
OBS=  15  at  0=  10  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 


Figure  27.  Comparison  of  the  Bias  of  (il  between  Different  numbers  of  targets  with 
OBS  =  35  at  0=10  in  the  GAM/YVEI  Regression  Model:  OBS  =  the 
number  of  observers 
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Table  16.  Tendencies  of  Mean  Biasand  Mean  Square  Error  for  /J1  at  O  =  15  in  the 
GAM/ WEI  Regression  Model 


M.B(M.S.E)  of  fa 

15TGT 

25  TGT 

35  TGT 

15  OBS 

.008(.021) 

-.Q06(.016) 

.0080011) 

25  OBS 

-.021  (.01 6) 

-.0040008) 

-.005(.006) 

35  OBS 

.015(.015) 

.008(.008) 

.001 (.004) 

Figure  28.  Comparison  of  the  Bias  of  /Jl  between  Different  numbers  of  targets  with 
OBS=  15  iX  0=  15  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers  t- 
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Figure  29.  Comparison  of  the  Bias  of  pi  between  Different  numbers  of  targets  with 
OBS  =  25  at  0=  15  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 


Figure  30.  Comparison  of  the  Bias  of  /ll  between  Different  numbers  of  targets  with 
OBS  =  35  at  0=  15  in  the  GAM/WEI  Regression  Model;  OBS  -  the 
number  of  observers 
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•4.  TABLES  AND  HISTOGRAMS  FOR  /?2 


Table  17;  Mean.  BiasjMean  Square  Error  and  Std.  Error  for  p2  at  O  =  10  in  the 
GAM/WEI  Regression  Model _ _ , 


Number  of 

*’  *  A 

Pi 

Observers 

Targets 

M.B(S.E) 

M.S.E(S.E) 

15 

,  .0167(.0097) 

,0096(;0012) 

15 

25 

,0057(;0075) 

.0056(.0008) 

35 

-,0097(.0053) 

.0029(.0004) 

15 

-.0051(10073) 

.0053(.0008) 

25 

25 

.0040(.0054) 

.0029(.0004) 

35 

,0020(.0046) 

,0021(.0003) 

15 

.0002(.0062) 

.0038(.0006) 

15 

25 

.008 3 (.0044) 

,0020(.0003) . 

35 

-.0046(.0033) 

.OOll(.OOOl) 

Table  18.  Mean  Bias, Mean  Square  Error  and  Std.  Error  for  P2  at  O  =  15  in  the 
_ GAM/WEI  Regression  Model _ _ 


Number  of 

Number  of 

A 

Pi 

Observers 

Targets 

M.B(S.E) 

M.S.E(S.E) 

15 

.0185(.0091) 

.008400012) 

15 

25 

.0Q55(.0070) 

.0049(.0006) 

35 

-.0124(.0052) 

.002S(.0004) 

15 

-.001500069) 

.004600006) 

25 

25 

.007200050) 

,0025(.0003) 

35 

.0024^.0042) 

,0018(.0002) 

15 

-.0016(.0057) 

.003200005) 

35 

25 

.0086(.0045) 

,0020(.0003) 

35 

-.0007(.0032) 

.OOIO(.OOOI) 
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Table  19.  Tendencies  of  Mean  Bias  and  Mean  Square  Error  for  {S2  at  0  =  10  in  the 


GAM/ WEI  Regression  Model 


I  M.B(M.S.E)  of  /h 

15  TGT 

25  TGT 

35  TGT 

15  OBS 

.017(.010) 

.0060006) 

-.0100003) 

25  OBS 

-.0050005) 

.0040003) 

.0020002) 

35  OBS 

.0000004) 

.0080002) 

-.0050001) 

OBS  =  15  ,  0  =  10  (  UC  :  76  PERCENT  ) 

N  -  15  N  *  25 


N  «  35 


Figure  31.  Comparison  of  the  Bias  of  p2  between  Different  numbers  of  targets  with 
OBS=  15  at  0=  10  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 
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Figure  32.  Comparison  of  the  Bias  of  pi  between  Different  numbers  of  targets  with 
OBS  =  25  at  0=  10  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 
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Figure  33.  Comparison  of  the  Bias  of  /J2  between  Different  numbers  of  targets  with 
OBS  =  35  at  0=  10  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 


Table  20.  Tendencies  of  Mean  Bias  and  Mean  Square  Error  for  /J2  at  O  =  15  in  the 


QAM/WEI  Regression  Model 


M.B(M.S.E)  of  & 

15  TGT 

25  TGT 

35  TGT 

15  OBS 

.018(.008) 

,006(.005) 

-.012(.003) 

25  OBS 

-.001(:005) 

.007(.002) 

.002(.002) 

35  OBS 

-.002(.003) 

.009(.002) 

-.OOl(.OOl) 

OBS.  =  15  ,  0  =  15  (  UC  :  85  PERCENT  ) 

N  -  15  N  -  25 


Figure  34.  Comparison  of  the  Bias  of  /?2  between  Different  numbers  of  targets  with 
OBS=  15  at  0=  15  in  the  GAM/WEI  Regression  Model:  013S  =  the 
number  of  observers 
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Figure  35.  Comparison  of  the  Bias  of  (12  between  Different  numbers  of  targets  with 
OBS  =  25  at  0=  15  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 


OBS  =  35  ,  0.  =  15  (  UG  :  85  PERCENT  ) 

N  =  15  N  «  25 


Figure  36.  Comparison  of  the  Bias  of  /?2  between  Different  numbers  of  targets  with 
OBS  =  35  at  0=  15  in  the  GAM/WEI  Regression  Model:  OBS  =  the 
number  of  observers 
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5.  TABLES  AND  HISTOGRAMS  FOR  $ 


Table  21.  Mean  Bias,  Mean  Square  Error  and  Std.  Error  for  £  at  O  =  10  in  the 
_ GAM/ WEI  Regression  Model _  _ 


Number  of 
Observers 

Number  of 
Targets 

A 

$ 

M.B(S.E) 

M.S.E(S.E) 

15 

15 

-.0318(.0067) 

.005400007) 

25 

-.0158(.0056) 

.003300005) 

35 

-.013100048) 

.002400003) 

25 

15 

-.022600053) 

.003300005) 

25 

-.007S00036) 

.001300002) 

35 

-.002500038) 

.001400002) 

35 

15 

-.019600044) 

.002300003) 

25 

-.006500036) 

.001300002) 

35 

-.007100028) 

.000800001) 

Table  22.  Mean  Bias,  Mean  Square  Error  and  Std.  Error  for  f  at  O  =  15  in  the 
GAM/WEI  Regression  Model  _ _ 


Number  of 
Observers 

Number  of 
Targets 

A 

M.B(S.E) 

M.S.E(S.E) 

15 

15 

-.034900064) 

.005200011) 

25 

-.014700051) 

.002700005) 

35 

-.012900041) 

.001800002) 

25 

15 

-.018100047) 

.002500003) 

25 

-.008900037) 

.001400002) 

35 

-.002700035) 

.001200002) 

35 

15 

-.019200041) 

.002000003) 

25 

-.004900033) 

.0011  (.0002) 

35 

-.00S300025) 

.0007(.0001) 
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Table  23.  Tendencies  of  Mean  Bias  and  Mean  Square  Error  for  l  at  O  =  10  in  the 


GAM /WEI  Regression  Model 


Figure  37.  Comparison  of  the  Bias  of  c  between  Different  numbers  of  Obse  rvers 
with  TGT=  15  at  0=  10  in  the  GAM/WEI  Regression  Model:  TGT 
=  the  number  of  tarsets. 

V1 


Figure  39.  Comparison  of  the  Bias  of  g  between  Different  numbers  of  Obse  rvers 
with  TGT=35  at  0=  10  in  the  GAM/WEI  Regression  Model:  TGI 
=  the  number  of  targets. 


Table  24.  Tendencies  of  Mean  Bias  and  Mean  Square  Error  for  £  at  O  =  15  in  the 
GAM/ WEI  Regression  Model 


M.B(M.S.E)  of  1 

15  TGT 

25  TGT 

35  TGT 

15  0BS 

-.035(.005) 

-.0150003) 

-.0130002) 

25  OBS 

-.0180002) 

-.009(.001) 

-.0030001) 

35  OBS 

-.0190002) 

-.0050001) 

.  -.0080001) 

Figure  40.  Comparison  of  the  Bias  of  £  between  Different  numbers  of  Observers 
with  TGT=  15  at  0=  15  in  the  GAM/WEI  Regression  Model:  TGT 
=  the  number  of  targets. 
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Figure  41.  Comparison  of  the  Bias  of  g  between  Different  numbers  of  Observers 
with  TGT=25  at  0=15  in  the  GAM/WEI  Regression  Model:  TGT 
=  the  number  of  targets. 
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Figure  42.  Comparison  of  the  Bias  of  $  between  Different  numbers  of  Observers 
with  TGT=  35  at  0=15  in  the  GAM/WEI  Regression  Model:  TGT 
=  the  number  of  targets. 
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APPENDIX  C.  SIMULATION  PROGRAM  FOR  THE  ESTIMATE  OF  tj 
IN  THE  SINGLE  PARAMETER  GAMMA  MODEL 


V  R*SIMULA1  iMiNiOiliJiKiLiETAi OUT; PAGE 5 A VUCLl MB ; SEMB 5 MSE ; SEMSE 

i  RRR 

Cl]  ARflflflflflflflfl  KEY  GLOBAL  VARIABES  :  BOOK 1  ;  BOOK 2  ;  UCL  HRflfiflflflflfi 

[2]  flflflflRflflRRflflflflflflflflflflflflflflflflflR  DATA  INPUT  RflflRRRRRRRRflRRfiflflflflflRflfl 

[3]  M*  15  25  35 

[4]  N*  15  25  35 

[5]  0 *  10  15 
C6]  ETA* 1 

[7]  I*J*K*L* 0 

C8]  RflflflRflRRflflflRflflRRflfl  START  LOOPING  THE  SIMULATION  1  rrrrrrrrrrrrr 

[9]  L00PI: 1*1+1 

[10]  LOOPJiJ*J+ 1 

[11]  L00PK:K*K+ 1 

[12]  RRflRflRRflflflflflflRRRflflR  RESET  THE  RANDOM  NUMBER  SEED  RflRRRRRRflflflRR 

[13]  []+0[I3tM[«7],JV[Ii:] 

[14]  Di?£<-466801743 

[15]  L00PLiL*L+ 1 

[16]  RflRRRRRflflRRflflRflfl  CALL  GAMMA  ESTIMATE  FUNCTION  RRRRRRRRRRRRflflRR 

[17]  RRR*QRL 

„  [18]  OUT*GAMMA(Oin  ,MIJ1  ,NIK1) 

[19]  rrparrr  RECORD  UCL  AND  ALL  THE  ESTIMATES  IN  REPLICATIONS  flflRfifl 

[20]  *UL=1)  ,L>1')/INTZ1,G01 

[21]  INTZliPAGE*(l,pOUT)p(OUT-0,ETA ) 

[22]  *LOOPL 

[23]  GOl: PA GE*PA GE , [1] (OUT-O , ETA ) 

[24]  TESTLl 

[25]  *(L<100)/LOOPL 

[26]  fl fl R fi r fl R  COMPUTE  AVE.  UN CENSORED  LEVEL  AND  ALL  THE  STATISTICS  Rfi 

[27]  AVUCL* ((  +  / PAGE l  •,  13  )+L) 

[28]  MB* (+ /PAGE [} 23 )+L 

[29]  SEMB*( (+/(PAGE[; 23 -MB )*2)+(LxL-2))*0. 5 

[30]  MSE*(+/PAGE[;23*2)+L 

[31]  SEMSE*((+/((PAGE[;23*2)-MSE)*2)+(LxL-2))*0.5 

[32]  L* 0 

[33]  RRRRRRRflfl  KEEP  ALL  THE  INFORMATIONS  FOR  GRAPHICS  AND  TABLES  flfl 

[34]  *{((,IxJxK)  =  l),  (IxJxK)*l)/INTZ2,G02 

[35]  INTZ2 : 

[36]  BOOKl* (l,pP AGE [; 2] )pPAGEi\23 

[37]  B00K2 * ( 1 , 4 ) p (MB, SEMB , MSE , SEMSE ) 

-  [38]  UCL*AVUCL 

[39]  +TESTK 

[40]  G02\ 

,  [41]  BOOK1*BOOK1, [1]  PAGE  Lx  23 

[42]  BOOI(2*BOOK2  ,  [  1  ]  (MB ,  SEMB  ,  MSE ,  SEMSE  ) 


91 


[43]  UCL+UCL ,AVUCL 

[44]  TESTK: 

[45]  HK<pN)/LOOPK 

[46] -  K+  0 

[47]  +U<pM)/LOOPj 

[48]  J+0 

[4.9]  +(I<pO)/LOOPI 
V 


V  R*GAMMA  INPUT',MiNiOiETAiETAOiXIiXIOiAiWH;MU;UiYiDELTAiCiXBiRR 
i Ml;  M2  ; T iDATA  \K%FL% FM.\FR ;LB;MP;RB;S 
Cl]  flflflflflftflflflftflfl  SIMULATION  INPUTS  AND  THEORETICAL  PARAMETERS  rrrr 
,  C2]  M+INPUT 121 

[3]  N+INPUT 131 
14]  O+XM ,N)p INPUT [ 1 ] 

C53  ETA* 1 

[6]  MU+n.2 

[7]  XlQ+§(N,M)pXI+MpO 

[8]  RflRRRRRflflflflflflRflflflR  RANDOM  NUMBER  GENERATIONS  rrrrrrrrrrrrrrrrr 
C 93  A+§(NtM)pM  GAMRAND ( (,*ETA ) ,  **ETA  ) 

CIO]  W+(M,N)p(MxN)  EXPRAND  1 

Cl 1]  flflflflfiflflflfiflflRR  DETECTION  TIME  AND  UNCENSORING  DATA  rrrrrrrrrrrrr 
C12]  U+MUxWW*A)*(*XI  0) 

Cl  3]  A+M+\  0 

ci4]  x+(®{/)L(®o) 

G 1 5]  0+/C2]  DELTA* (®£/)<;(®0) 

[16]  0* l0 

Cl 7]  XB+toMU 

CIS]  RRRRflRRflflflRR  INITIAL  CONDITION  FOR  GAMMA  PARAMETER  tETA  RRRflflR 
C 1 9  ]  /?iMY-XB)x*ClxXlO) 

C  2  0 ]  Ml«-“0 .5772 

C21]  M2*(,+/+/DELTAxRR*2)*(+/C) 

C22]  rao<-o 

C/23]  +(0>T+(  (  (ol  )*2  )i~6  )+(M2*2  )-  (Ml*2  )  )/BISECTION 

J  [24]  ETA0+~lx®T 

[25]  RRflRR  FIND  BOUNDS  AND  BISECTION  SEARCH  FOR  GAMMA  PARAMETER , ETA  fifl 
C  2  6 ]  BISECTION lK+0 
*  C  2  7 ]  DATA+(S*+/121*RR),C 

C 2 8 3  FM+ETAO  FVALUE  DATA 

C 2 9 ]  RRRARRRRflRflflflfi  FIND  THE  RIGHT  AND  LEFT  BOUND  FOR  ETA  rrrrrrrrrrr 

[30]  BOUND iK+K+1 

[31]  FL+(.LB+ETA0-KX0.5)  FVALUE  DATA 

[32]  FR*(RB*ETAO+KxO , 5 )  FVALUE  DATA 

[33]  *(((xFM)x(xFL))<0)/SETRB 

C34]  •*■(  ( (xFM)x (xFR) )<0 )/SETLB 

[35]  +BOUND 

[36]  SETRB :RB*ETA0 
[373  ->BISECT 

[38]  SETLB : LB+ETA 0 

[39]  rrrrarrrrrrrrrr  PERFORM  THE  BISECTION  SEARCH  FOR  ETA  aararaaaaa 

[40]  BISECT i 

[41]  FM+(MP+(LB+RB)*2)  FVALUE  DATA 

[42]  FL*LB  FVALUE  DATA 

[43]  +URB-LB)<(  \  Q.01*MP))/END 

1  [44]  *((((xFM)x(xFL))<0), (((xFM)x(xFL))>0))/SNAPR,SWAPL 

[45]  SNAFR : RB+MP 

[46]  +BISECT 

■  [47]  SNAPL : LB<-MP 

[48]  ^BISECT 
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[49]  AAAAAAAAAAAAAAAAAAAAAAA  OUTPUT  THE  ESTIMATES  AAAAAAAAAAAAAAAAA 

[50]  END l 

[51]  ETA0+(RB+LB)*2 

[5  2  ]  jM(+/C)  )  ,ETA0 

V 


V  R+ETA  FVALUE  DATA\M\S\C\SUBF\I 

[1] 

[2]  S+M+DATA 

[3]  WtMM 

[4]  SUBF+O 

[5]  i+o 

[6]  L0OPM+Z+1 

[7]  *(ClTl=0)/CHECKI 

[8]  SUBF+SUBF+(+/*(*ETA)+(\CU)  )-l) 

[9]  CHECKI : 

[10]  ^ (I<M) /LOOP 

[u]  Pf-c*m)xs£/BF+(+/m+<(s-c)^s+*EW)+“ix®(s+*m)) 

V 


APPENDIX  D.  SIMULATION  PROGRAM  FOR  THE  ESTIMATES  OF 
ALL  PARAMETERS  IN  GAMMA/WEIBULL  REGRESSION  MODEL 


C-l] 

123 

C3] 

[4] 

C'5] 

C63 

C7] 

[8] 

X93 

CIO] 

Cll] 

C12] 

C13] 

Cl  4] 

Ci;5] 

C16] 

C17] 

C18] 

C19] 

[20] 

C21] 

C  2  2  ] 

C  2  3  ] 

C24] 

C25] 

C  26  ] 

C  2  7  ] 

C28] 

C29] 

C30] 

C31] 

C32] 

C33] 

C34] 

C35] 

C36] 

C37] 

C38] 

C  3  9  3 

C40] 

C41] 

C42] 


V  R+SIMULA2iMiNiOiIiJiKiLiXiEiETAiBOiBUB2;XhOUTiPAGE;AVUCEXiMB 
iSEMBiMSEiSEMSEiRRR 

flflflflflflfiflflflfl  KEY  GLOBAL  VARIABES  i  B00K1  ;  B00K2  ;  UCLEX  RflnfiRAfififl 
RRRRRRRRRRRRRRRRRRRRRRRRRRR  DATA  INPUT  RRRRRRRRRRRRRRRRRRRRRR 


15  25  35 
N+  15  25  35 
0+  10  15 
ETA+1 


B 0+0,8 

2?l«-“0 » 2 


B2«-0 . 5 
XI«-0 


RflRflflRRflflflflRRRflfl  START  LOOPING  THE  SIMULATION  ARfiflfiflflflARRflfiAflRR 
I +J  <-K+L+E+X+  0 
LOOPI :I+I+1 
L00PJ'.J+J+1 
LOOPKiK+K+1 

RflflflflRRflRflRRRRRRRR  RESET  THE  RANDOM  NUMBER  SEED  RRRfiRfiRRRRfifiRRA 

□+0CI],MC«7],tfCX] 

□££<•466801743 

L00PL',L+L+1 

RRRflflRflflRRRR  CALL  GAMMA/WEIBULL  REGRESSION  FUNCTION  RflflflflflflflRRR 
RRR+URL 

OUT+GAMWEI(.Oin  ,M£J1  ,NIK1 ) 

RRfififififlflflfiR  SKIP  THE  INSUFFICIENT  VARIABILITY  CASE  FOR  GAMMA  fifl 
+  (.OUT C 2 ]  *“ 9 9 9 9 9  ) /CHECKETA 
E+E+ 1 


L+L- 1 


+TESTL 


rrrrrr  RECORD  THE  RANDOM  SEED  FOR  THE  ILL  CONDITION  OF  ETA  flflfl 
CHECKETA : 

(002^2]  *“100000  )/GOGO 
D+'BAD  ETA » 

□ +RRR 
L+L-l 


+TESTL 


RflRfiRRRRfl  RECORD  UCL  AND  ALL  THE  ESTIMATES  IN  REPLICATIONS  RRfl 
GO  GO  • 

♦( (£=1 ) ,L>1 )/INTZl ,G01 

INTZ1 1 PAGE* (1 ,  (,pOUT)~l  )p  (  (~1$OUT)-0 , ETA, B0  ,B1  ,B2  ,XI) 

X+X+C 1+OUT) 

+LOOPL 

G01 : PAGE+PAGE ,  Cl ] ( ( “1 *OUT ) - 0 , ETA , B0 , B1 , B2 , XI ) 

X+X+CHOUT) 
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[43]  RRRflflflflfl  RECORD  THE  RANDOM  SEED  FOR  THE  ILL  CASE  OF  XI  RflfifiRRfifl 

[44]  +tCl*OUT)-Q)ATESFL. 

[4.5]  0  +'BAD  XI ' 

[46] 

[47] ,  flflflflflfiflflflflRflflflflfiflflRRflRflflflflflflflRPRRRRRflflflflPRRSRfiRRPRRfiRfi  RR  RRRRflflPRfl 

[48]  TEST Li 

[49]  +(L<100)ALOOPL 

[50]  flflfiflflflfl  COMPUTE  AVE,  UN CENSOR*  >  LEVEL  AND  ALL  THE  STATISTICS  fl 

[51]  AVUCEX* ((+/ PAGE l ; 1] )*L),E,X 

[52]  MB*(+/l 1]  PAGE*  0  1  *PAGE)iL 

[53]  SEMB*((+/ 11]  (PAGE-  (pPAGE)pMB  )*2  )  +  ( L*L-2  ) )*0 . 5 

[54]  MSE«-(+/[l]  PAGE*2)*L 

[55]  SEMSE*(  (+/[1]  ( (PAGE*!  )-  (pPAGE)pMSE)*2  )*(LxL-2  )  )*0 . 5 

[56]  L+E+X+ 0 

[57!  RRflRRRflfi  KEEP  ALL  THE  INFORMATIONS  FOR  GRAPHICS  AND  TABLES  Rfifi 

[58]  +( ( (IxJxK)~ 1 ) , (IxJ*K)*l )/INTZ2,G02 

[59]  INTZ2 : 

[60]  BOOKl+(l ,pPAGE)pPAGE 

[61]  BOOK2+ ( 1 , 4 , pMB ) p (MB , SEMB , MSE , SEMSE ) 

[62]  UCLEX+  1  3  p AVUCEX 

[63]  +TESTK 

[64]  GO 2: 

[65]  BOOKl+BOOKl,ll]  PAGE 

[66]  BOOK2*BOOK2,Ll]  ((V  ,pMB)pMB ,  SEMB , MSE ,  SEMSE) 

[67]  UCLEX+UCLEX,l  1]  AVUCEX 

[68]  TESTK : 

[69]  +(K<pN)/LOOPK 

[70]  K* 0 

[71]  +(J<pM)/LOOPJ 

[72]  «7«-0 

[73]  *(I<pO)/LOOPI 
V 
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VGAMWEI  [D3V 

V  R*GAMWEI  INPUT iM;NiO; ETA ; ETAO  ;  B ;  £0 ;  B1 ;  £2 ;  XJ ;  XIO ;  ,5 ;  JiW;  Xl }  X2 
; MU •, Ui Ti DELTA ; C ;  XB ;RR;Wi'U 0  ; U t; U2;UU;Z;I \VAR ; Ml \M2\T\ DATA ; K 
iFL;FLO.;FM;FRiFRO;LB;MPiRBiCSiH;DH;DIFFXIiBKiSiJ;XiXXiBADXI 
111  flflflflflpftflflflfl  SIMULATION  INPUTS  AND  THEORETICAL  PARAMETERS  Rflfiflflflfl 
[23  M+INPUT121 

[3]  N+INPUTC31 

[4]  0+{M,N)p INPUT  111 

[5]  ETA+1 

t6]  B0<-0.8 

[73  £l<-~0.,2 

[83  B2+0.5 

[93  XIO*^(,N,M)pXI*MpO 

[103  fl.flRRRflRRflflflflflflflflfl  RANDOM  NUMBER  GENERATIONS  RRRflRRRflflflflflRRfififlflflfl 
[113  A+§(N,M)pM  GAMRAND((*ETA) ,**ETA) 

[123  WW+(M,N)p(MxN)  EXPRAND  1 
'C133  XI*  07,  AO  P  (M*N)  NORRAND  1  0.5 
[143  X2*(M,AOp(MxAO  NORRAND  2  1 

[153  RRflflflRflflRflflfl  DETECTION  TIME  AND  UNCENSORING  DATA  RflfiflflRRRRRRRRRR 
[163  MU+*  (£0+ (BlxXl  )+B2xX2  ) 

[173  U<rMU*(W*A)*(*XI  0) 

[133  A*NN+MU+\0 
[193  I*(©£OL (®a) 

[203  C++/L21  DELTA* (®U)4(®0) 

[213  0+\ 0 

[223  RRflflflRflRPRfl  INITIAL  ESTIMATES  WITHOUT  GAMMA  R.V.  INVOLVED  Rflflflflfl 

[233  XB*®U 

[243  1+ 0 

[253  INTZ : 1*1 +1 

[263  W<-(*  (Y-XB)x*  C”lxXI0  )  )*0 . 5 

[273  UO+((MtN)pl)xWx*Cl*XIO) 

[283  i/l^-Xlx^x*("lxXJ0) 

[293  U2+X2*Wx*Cl*XI0) 

[303  UU*§(3 ,MxN)p ( ,U0) , ( ,Ul) , ( ,U2) 

[313  U0*Ul*U2*\0 

[323  Z*§,  (WxXBxi<rixXI0))  +  (,rixDELTA')+W'K2)*W 

[333  N+\ 0 

[343  B*ZWU 

[353  Z*UU*\ 0 

[363  BK*B 

[373  XB*£[13+(B[23xXl)+fl[33xX2 

[383  Bi?-(y-XB)x*rixXJO) 

[393  XJO*§ (N ,M)pXI*XI+DIFFXI*(C-+/RRx  (  Cl*DELTA  )+*££))* 

( (~lxC)-+/(i?i?*2  )**RR) 

[403  *U<2)/INTZ 

[413  RflflfiRRRRflflRRR  VARIABILITY  CHECK  FOR  UNCENSORED  TIMES  RflRRRflfiflflfifl 
[423  VAR*(.+/+/Cl  +  (U*(*XB)')*(*~lxXI0))*2)H~l+MxN) 

[433  U*i 0 
[443  £ZVi0*~99999 

[453  X+Mpl 
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[46]  J+ 0 

[47]  +(VAR<1)/END 

[483-  fl  fl’fl  fi  fl.fl  fl  fl  fl  fl  fl  fl  ft  INITIAL  CONDITION  FOR  GAMMA  PARAMETER  ,ETA  flflfiflflflfl 

[49]  RR+a-XB)x*  ("IxXlO ) 

[50]  Ml<-“0.5772 

[513=  m<r(+/+/DELTA*RR*2)*(+/C) 

[523  ETAQ+O 

[533  +  (0>2’*(((ol)*2)*”6)+(M2*2MMl*2))/flE;iDy 

[543  ErAO-^lx®^ 

[553  READHBADXI+ 0 

[563  Rfl  LARGE  RECURSIVE  LOOP  INCLUDING  ESTIMATE  OF  GAMMA  PARAMETER , ETA 
[573  RECURSIVE iK+ 0 
[533  c7«-«7+l 

[593  DATA+(,S<r+'/L2l*RR),C 
[603  FM+ETAO  FVALUE  DATA 

[613  AflflRRflfl  BISECTION  SEARCH  METHOD  APPLIED  TO  ESTIMATE  ETA  rrrrrrrr 
[623  BOUND XK+K+ 1 

[63]  FL+(LB+ETAO-KxO.$)  FVALUE  DATA 
[643  FR+(RB+ETAO+KxO.S) FVALUE  DATA 

[653  RflRRRRRRflRflRflflflflflR  DETECT  ILL  CONDITION  FOR  ETA  RRRRfiflfiflfiAflflARRA 
[663  +(K£10) /CHECKLRB 

[67]  -►('v(  ( (  |  FRO-FR)*Q  ,  5  )<0 . 01  )A (  (  |Fi?)<0 . 1 )) /CHECKLRB 

[68]  EM0«-~100000 

[69]  +END 

[70]  flflRRRflRRRRflRflRRflflflhRRRflRRflflRflflRflRRRRflRRflRflRRflRRRRRRRRRflRRflRRRRRR 

[71]  CHECKLRB X 

[72]  •>(  (  (xfW)x  (xFL)  )<0  )/SETRB 

[73]  •*(  (  (xFM)x  (xfTJ )  )<0  ) /SETLB 

[74]  FLO+FL 

[75]  FRO+FR 

[76]  ABOUND 

[77]  SETRB x RB<-ETA 0 

[78]  +BISECT 

[79]  SETLB xLB+ETAO 

[80]  BISECT x 

[81]  FM+(MP+(LB+RB)*2)  FVALUE  DATA 

[82]  FL+LB  FVALUE  DATA 

[83]  *»■(  (RB”LB  )<  (  i  0 .  OlxMP)  ) /REGRESS 

[84]  *(  (  (  (xFM)x  (xFL  )  )<0  )  ,  <  (  (xfM)x  (xFL)  )>0  )  )  /  SWAPR  ,SWAPL 

[85]  SWAPR XRB+HP 

[86]  +BISECT 

[87]  SNAPLxLB+MP 

[88]  +BISECT 

[89]  RRflfifl  REGRESSION  PROCEDURE  FOR  WEIBULL  SCALE  PARAMETERS , BETA ' S  RR 

[90]  REGRESS  X 

[91]  ETAQ<-(RB+LB)*2 

[92]  W+((CS+§(N,M)  p  (  ( C+*ETAO  )*S+*ETAO  )  )x*/?H)*0 . 5 

[93]  £/Ot-((M>W)pl)x^x*(”lxXIO) 

[94]  Ul+XlxNx*  (*"lxXJO  ) 

[95]  i/2<-X2xp/x*  (“IxXJO  ) 

[96]  UU+§( 3 ,MxN)p ( ,UQ ) , ( ,U1 ) , ( ,U2 ) 
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[97]  U.0*Ul+U2+\0 

•098]  Z+§  ,<W*XBx*  ClxXIO  )  )+  (  ClxDELTA  )+W* 2  )*W 

099]  B+muu 

01 00]  UU+Z+ lO 

C101]  Rfl-flflflflflflfl  NEWTON  PROCEDURE  FOR  WEIBULL  SHAPE  PARAMETER , X I « S  RRR 
-01.02]  XB+S[l]+(BC2]xXl)+BC3]xX2 
0103]  RR+(I-XB)x*Cl*XlO) 

0104]  B-<-(“lxC)  +  (+/02]  RRx(.ClxDELTA)+CSx*RR)') 

Ci05]  M*OxC)'-(+/C2]  CSx(*RR)*RR* 2) 

0106]  flflflflflflflflflflflflflflflR  DETECT  THE  ILL  CONDITION  OF  XI'S  RflflflflRflflRflflflfl 
0107]  XX+~(  (|B)<0.1)a((|BB)<0.01) 

0108]  +  ( ~ ( ( x /XX )  =  0 ) aB4BXI= 0 )  / NEWTONXI 

0109]  BADXIi+1 
CllO]  X+XX 
[111]  NEWTONXI i 

0112]  XI 0^  </v\  M  )  pXI+XI+D  JFFXIx-XXx  ( “  1  xB ) +BB 

0113]  rr  TEST  FOR  STOPPING  CRITERIA  AND  UPDATE  THE  VALUE  OF  VARIABLES 
0114]  +((x/( | (DIFFXI+XI) , (BK-B)*B )<0 . 01 )v (JS50 ) )/END 

0115]  BK+B 

0116]  BB-t-a-XB)x*(-lxXIO) 

0117],  •> RECURSIVE 

0118]  RRRRflflRRflfl  OUTPUT  UN CENSORED  LEVEL  AND  ALL  ESTIMATES  RflflRRRRRRfl 
0119]  END: 

0120]  B<-((+/C)*M),BM0,BCl],B[2],B[3],((+/XrxX)*  (+/X)),+/~X 

V 
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